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ABSTRACT 

We report on the X-ray evolution over the last as 9 Gyr of cosmic history (i.e., since z = 1 .4) of late-type galaxy 
populations in the Chandra Deep Field-North and Extended Chandra Deep Field-South (CDF-N and E-CDF-S, 
respectively; jointly CDFs) survey fields. Our late-type galaxy sample consists of 2568 galaxies, which were 
identified using rest-frame optical colors and HST morphologies. We utilized X-ray stacking analyses to in- 
vestigate the X-ray emission from these galaxies, emphasizing the contributions from normal galaxies that are 
not dominated by active galactic nuclei (AGNs). Over this redshift range, we find significant increases (factors 
of «5-10) in the X-ray-to-optical mean luminosity ratio (LxILb) and the X-ray-to-stellar-mass mean ratio 
(LxIM*) for galaxy populations selected by Lb and M*, respectively. When analyzing galaxy samples selected 
via SFR, we find that the mean X-ray-to-SFR ratio (Lx/SFR) is consistent with being constant over the entire 
redshift range for galaxies with SFR = 1-100 M© yr" 1 , thus demonstrating that X-ray emission can be used as 
a robust indicator of star-formation activity out to z ~ 1 .4. We find that the star-formation activity (as traced by 
X-ray luminosity) per unit stellar mass in a given redshift bin increases with decreasing stellar mass over the 
redshift range z = 0.2-1, which is consistent with previous studies of how star-formation activity depends on 
stellar mass. Finally, we extend our X-ray analyses to Lyman break galaxies at z ~ 3 and estimate that LyjLs 
at z ~ 3 is similar to its value at z = 1 .4. 

Subject headings: cosmology: observations — surveys — galaxies: normal — galaxies: spiral — galaxies: 
star formation galaxies: active — galaxies: — X-rays: galaxies — X-rays: general 



1. INTRODUCTION 

Investigations focusing on global changes in star-formation 
activity and stellar-mass build-up in field galaxies have pro- 
vided significant insight into the physical evolution of galax- 
ies and their constituent stellar populations. It has now been 
well-established that the global star-formation rate density has 
declined by roughly an order of magnitude since z ~ 1—1.5 
(e.g., Lilly etal. 1996; Madau etal. 1996; Steidel etal. 1999; 
Hopkins etal. 2004; Perez-Gonzalez etal. 2005; Schimi- 
novich etal. 2005; Colbert etal. 2006). Recent investiga- 
tions into the details of this evolution have shown that the 
star-formation history of a given galaxy depends strongly on 
its stellar mass (e.g., Cowie etal. 1996; Juneau etal. 2005; 
Bundy etal. 2006; Noeske etal. 2007 a,b; Zheng etal. 2007); 
the peak star-formation epoch for the most massive galaxies 
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occured earlier in cosmic history than it did for galaxies with 
lower masses. 

X-ray studies of normal late-type galaxies (i.e., those 
that are not dominated by luminous active galactic nuclei 
[AGNs]) have shown that X-ray emission provides a useful, 
relatively-unobscured measure of star-formation activity (e.g., 
Bauer etal. 2002a; Cohen 2003; Ranalli etal. 2003; Colbert 
etal. 2004; Grimm etal. 2003; Gilfanov etal. 2004a; Per- 
sic etal. 2004; Persic & Rephaeli 2007; however, see Barger 
etal. 2007). In normal galaxies, X-ray emission originates 
from X-ray binaries, supernovae, supernova remnants, hot 
(wO.2-1 keV) interstellar gas, and O-stars (see, e.g., Fab- 
biano 1989, 2006 for reviews). Sensitive Chandra and XMM- 
Newton studies of local late-type galaxies have shown that 
high-mass X-ray binaries (HMXBs) and low-mass X-ray bi- 
naries (LMXBs) typically dominate the total non-nuclear 
X-ray power output (e.g., Zezas et al. 2002; Bauer et al. 2003; 
Soria & Wu 2003; Swartz et al. 2003; Jenkins et al. 2005; Kil- 
gard etal. 2005; however, see, e.g., Doane etal. 2004). Ob- 
servations indicate that the integrated X-ray emission from 
HMXB and LMXB populations trace galaxy star-formation 
rate (SFR) and stellar mass (M*), respectively. For example, 
using Chandra observations of 32 local galaxies, Colbert et al. 
(2004) found that the summed 0.3-8 keV non-nuclear point- 
source emission from a given galaxy (Lxp) can be approxi- 
mated as Lxp ~ a + (3 SFR, where a and (3 are constants. 
Therefore, galaxies having relatively high star-formation rates 
per unit mass (specific star-formation rates, SSFRs) gener- 
ally have dominant X-ray point-source contributions from 
HMXBs (e.g., late-type star-forming galaxies), while those 
with relatively low SSFRs have point-source emission primar- 
ily from LMXBs (e.g., massive early-type galaxies). 

If the X-ray binary populations are similarly dominating the 
normal-galaxy X-ray power output over a significant fraction 
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of cosmic time, then there should be a rapid increase in the 
globally-averaged X-ray luminosity of normal star-forming 
galaxies with cosmic look-back time in response to the in- 
creasing global star-formation rate density (e.g., Ghosh & 
White 2001). In this scenario, HMXBs trace the immedi- 
ate star-formation rate of a galaxy and LMXBs trace its star- 
formation history with a lag of a few Gyr. With the advent of 
deep Chandra and XMM-Newton surveys (see, e.g., Brandt 
& Hasinger 2005 for a review), it has become possible to 
study the X-ray properties of normal galaxies out to z 1 
andz w 0.3, respectively (see, e.g., Hornschemeier et al. 2000, 
2003; Alexander etal. 2002; Georgakakis etal. 2007; Geor- 
gantopoulos et al. 2005; Lehmer et al. 2006, 2007; Kim et al. 
2006; Tzanavaris etal. 2006; Rosa-Gonzalez etal. 2007). Ev- 
idence for a global increase in the X-ray activity with red- 
shift for normal late-type galaxies has since been mounting. 
Studies of the X-ray luminosity functions of X-ray-detected 
normal galaxies have found that evolves as (1 + z) 15 " 3 
over the redshift range z ~ 0-1.4 (e.g., Norman etal. 2004; 
Georgakakis et al. 2006; Ptak et al. 2007; Tzanavaris & Geor- 
gantopolous 2008). Additionally, X-ray stacking analyses 
have enabled investigations of more representative optically- 
selected galaxy populations over the majority of cosmic his- 
tory (z ~ 0.1—4; e.g., Brandt etal. 2001; Hornschemeier etal. 
2002; Nandra etal. 2002; Georgakakis etal. 2003; Reddy & 
Steidel 2004; Laird etal. 2005, 2006; Lehmer etal. 2005a). 
These studies have found that the average X-ray luminosi- 
ties of normal late-type galaxies increases with redshift out 
to z ~ 1.4-3. For example, using a «1 Ms exposure of a sub- 
region within the Chandra Deep Field-North, Hornschemeier 
etal. (2002, hereafter H02) tentatively observed a factor of 
r*2-3 increase in Lx/L B from z = to 1.4 for L* B galaxies. 
Despite these promising initial constraints, the details of the 
X-ray evolution of normal late-type galaxy populations, in- 
cluding dependences on the physical properties of galaxies 
(e.g., optical luminosity, stellar mass, environment, and star- 
formation rate), have remained unexplored. 

In this paper, we aim to improve significantly upon con- 
straints for the X-ray evolution of normal late-type galax- 
ies (e.g., H02). We study for the first time how the X-ray 
properties of late-type field galaxies evolve as a function 
of optical luminosity, stellar mass, and star-formation rate 
over the redshift range of z = 0-1.4. We construct late-type 
galaxy samples located in two of the most well-studied ex- 
tragalactic X-ray survey fields, the f«2 Ms Chandra Deep 
Field-North (CDF-N; Alexander et al. 2003) and the Extended 
Chandra Deep Field-South (E-CDF-S), which is composed 
of the central wl Ms Chandra Deep Field-South (CDF-S; 
Giacconi etal. 2002) and four flanking w250 ks Chandra 
observations (Lehmer etal. 2005b). These Chandra Deep 
Fields (hereafter CDFs) reach 0.5-2 keV detection limits of 
w2.5 x 10~ 17 ergs cm" 2 s" 1 in the most sensitive regions and 
^3 x 10~ 16 ergs cm" 2 s" 1 over the majority of the CDFs; 
these levels are sufficient to detect moderately-powerful X-ray 
sources (Lo.5-2 keV ^ 10 41 ' 5 ergs s" 1 ) at z = 1 .4 and z = 0.6, re- 
spectively. Therefore, the CDFs comprise an unprecedented 
data set for effectively studying the X-ray emission and evolu- 
tion of cosmologically distant normal galaxies with minimal 
contamination from powerful AGNs. 

The Galactic column densities are 1.3 x 10 20 cm" 2 for 
the CDF-N (Lockman 2004) and 8.8 x 10 19 cm" 2 for the 
E-CDF-S (Stark etal. 1992). All of the X-ray fluxes and lu- 
minosities quoted throughout this paper have been corrected 
for Galactic absorption. Unless stated otherwise, we quote 



optical magnitudes based upon the AB magnitude system 
for HST photometry and the Vega magnitude system every- 
where else. In the X-ray band, we make use of three standard 
bandpasses: 0.5-2 keV (soft band [SB]), 2-8 keV (hard band 
[HB]), 0.5-8 keV (full band [FB]). Throughout this paper, we 
make estimates of stellar mass and star-formation rates us- 
ing a Kroupa (2001) initial mass function (IMF); when mak- 
ing comparisons between these estimates and those quoted in 
other studies, we have adjusted all values to correspond to our 
adopted IMF. H = 70 km s" 1 Mpc" 1 , fl M = 0.3, and fl A = 0.7 
are adopted throughout this paper (e.g., Spergel etal. 2003), 
which imply a look-back time of 7.7 Gyr at z = 1 . 

2. LATE-TYPE GALAXY SAMPLE SELECTION 

We constructed an optically selected sample of late-type 
galaxies within the CDFs to use for our subsequent anal- 
yses. We restricted our galaxy selection to regions of the 
CDFs where Hubble Space Telescope (HST) observations 
were available to allow the best possible morphological clas- 
sifications. The HST observations in the CDFs have been 
carried out via the Great Observatories Origins Deep Sur- 
vey (GOODS; Giavalisco etal. 2004a) and Galaxy Evolution 
from Morphology and SEDs (GEMS; Rix etal. 2004; Cald- 
well etal. 2005) programs; these surveys cover «90% of the 
Chandra-observed regions of the CDFs with the Advanced 
Camera for Surveys (ACS). In the GOODS and GEMS re- 
gions, photometry was available for four (B435, V606, '775, and 
Z850) and two (Veoe and z&so) ACS passbands, respectively. 

2.1. Galaxy Selection Footprint 

We began building our sample by selecting all galaxies hav- 
ing Z850 < 23. This initial selection criterion was motivated by 
(1) the availability of deep z 850 band imaging over all of the 
CDFs, (2) the fact that the Z&50 emission probes mass-tracing 
rest-frame optical light redward of the 4000 A break for galax- 
ies at z 1, which constitutes the large majority of our sam- 
ple, (3) the availability of reliable spectroscopic and photo- 
metric redshifts (see details below), and (4) the high multi- 
wavelength detection fractions for these sources, which allow 
us to determine informative rest-frame spectral energy distri- 
butions (SEDs). In order to isolate most effectively distant 
X-ray-detected AGNs, we further restricted our sample to in- 
clude only galaxies that were located in the most sensitive 
areas of the CDFs where the Chandra point-spread function 
(PSF) was small. We therefore chose to include sources hav- 
ing optical positions that were within 6.'0 of at least one of 
the six Chandra aimpoints in the CDFs 12 ; the corresponding 
total areal footprint is «0.18 deg 2 . Futhermore, we removed 
obvious Galactic stars that were identified via optical spectral 
properties or the presence of obvious diffraction spikes in the 
Z850 band images. Under these restrictions, we found 6905 
galaxies. 

2.2. Redshifts 

We cross-correlated our initial sample of 6905 galaxies 
with available spectroscopic and photometric redshift cata- 
logs (e.g., Barger etal. 2003; Le Fevre etal. 2004; Szokoly 
etal. 2004; Wirth etal. 2004; Wolf etal. 2004; Mobasher 
etal. 2004; Mignoli etal. 2005; Vanzella etal. 2005, 2006; 
Grazian et al. 2006; Ravikumar et al. 2007; Silverman et al. 

12 For the CDF aimpoints, see Tables 1 and Al of Alexander et al. (2003) 
for the ^2 Ms CDF-N and ss 1 Ms CDF-S, respectively, as well as Table 1 of 
Lehmer et al. (2005b) for the «250 ks E-CDF-S. 




FIG. 1. — Rest-frame U — V color versus absolute V-band magnitude My (color-magnitude diagrams) for the 3259 z ~ 0-1.4 galaxies with Z850 < 23 that 
were within 6f0 of at least one CDF aimpoint. Each panel shows the color-magnitude relation for a given redshift bin (annotated in the upper left-hand corners). 
The dashed line in each panel represents the estimated division between red (large U - V values) and blue (small U — V values) galaxy populations, which was 
estimated using equation 1 and the median redshift of the galaxies in each bin. Open symbols highlight red-sequence and blue-cloud galaxies that were visually 
reclassified as late-type (open circles) and early-type (open squares) galaxies, respectively. The inset histogram in each panel shows the distribution of rest-frame 
U—V. The vertical dotted line in each inset plot indicates the estimated division between red and blue galaxy populations, which was calculated using equation 1 
and the median redshift and My for galaxies in each redshift bin (see § 2.3). 



2008). All galaxies that did not have spectroscopic redshifts 
were located in the E-CDF-S where highly accurate (median 
Sz/l+z ~ 0.02 for galaxies with z$so < 23) photometric red- 
shifts were available via COMBO-17 (Classifying Objects by 
Medium-Band Observations in 17 Filters; Wolf etal. 2004). 
In total, 6683 («97%) of our sources had either spectroscopic 
or photometric redshifts. Visual inspection of the 222 galaxies 
in the E-CDF-S without redshifts indicated that these sources 
were mainly faint galaxies near bright stars, as well as a hand- 
ful of sources that were subgalactic features within relatively 
nearby galaxies. 

Whenever possible, we adopted spectroscopic redshifts as 
the most accurate redshifts for our galaxies. Using the red- 
shift information, we filtered our sample to include only 
sources with z < 1.4 in the w2 Ms CDF-N, z < 1 in the 
«1 Ms CDF-S, and z < 0.6 in the «250 ks regions of the 
E-CDF-S; these redshift limits represent the largest distances 
at which we would expect to identify moderately luminous 
(Lq.5-2 kev ~ 10 41 ' 5 ergs s" 1 ) AGNs effectively in each respec- 



tive field. In total, 3259 galaxies remained after filtering our 
sample based on redshift properties. We used spectroscopic 
redshifts for 1351 (w41%) galaxies and photometric redshifts 
for the remaining 1908 («59%) galaxies. 

2.3. Rest-Frame Color and Morphological Selection 

The optical-color distribution for field galaxies has been 
shown to be bimodal, separating "red" and "blue" galaxy pop- 
ulations (e.g., Strateva etal. 2001; Hogg etal. 2002a; Blan- 
ton etal. 2003; Baldry etal. 2004). Studies of the color- 
magnitude relation for distant galaxy populations have shown 
that this color bimodality is observed to persist out to at least 
z w 1-1 .5 (see, e.g., Bell et al. 2004a; Faber et al. 2007; Labbe 
et al. 2007), thus providing an excellent quantifiable means for 
separating late-type and early-type galaxy populations. We 
therefore filtered our galaxy sample to include only sources 
that had blue rest-frame optical colors, as expected for late- 
type galaxies that contain young stellar populations. To es- 
timate rest-frame U — V colors and absolute U, B, and V 
band magnitudes (M v , M B , and M v , respectively) for each 
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of the 3259 galaxies in our sample, we utilized the available 
optical/near-IR data. 

For the CDF-N, we used the U, B, V, R, I, z' , and HK' pho- 
tometric catalogs from Capak et al. (2004), as well as Spitzer 
IRAC imaging (3.6, 4.5, 5.8, and 8.0^m; Fazio etal. 2004) 
from the GOODS (Dickinson etal., in preparation). For the 
E-CDF-S, we used (1) the 17-bandpass photometry available 
through COMBO- 17, (2) J and K s imaging from the ESO 
Imaging Survey (Olsen etal. 2006), (3) J, H, and K from 
the MUSYC collaboration (Taylor etal. in preparation; see 
also, Moy etal. 2003; Gawiser etal. 2006), and (4) Spitzer 
IRAC data from the GOODS and SIMPLE teams (Dickin- 
son etal., in preparation; van Dokkum etal., in preparation). 
Using these data, we constructed a rest-frame near-UV-to- 
near-IR SED for each galaxy. For these SEDs, the rest-frame 
U, B, and V filters are well-sampled at all relevant redshifts, 
with the exception of sources at z rt 0.7 where the wavelength 
range of the available data is sparse at rest-frame V band. For 
sources at z ^ 0.7, we linearly interpolated our SED to cover 
the V filter. We convolved these SEDs with Johnson U, B, and 
V filter curves and computed rest-frame absolute magnitudes 
for each respective filter following equation 5 of Hogg et al. 
(2002b). For sources in the E-CDF-S, these computed abso- 
lute magnitudes are consistent with those presented by Wolf 
etal. (2004). 

In Figure 1 (small filled circles), we present rest-frame 
U—V colors versus My for our galaxies in six redshift ranges. 
For clarity, we also show inset histograms giving the distri- 
bution of rest-frame U — V colors for each redshift interval. 
We utilized the rest-frame U — V color to divide roughly pop- 
ulations of red and blue galaxies. Use of the rest-frame U — V 
color was motivated by Bell et al. (2004a), who note that the U 
and V bandpass pair straddle the 4000 A break, which is par- 
ticularly sensitive to age and metallicity variations of galac- 
tic stellar populations. The dashed lines in Figure 1 show 
the empirically-determined redshift-dependent color divisions 
that separate blue and red galaxy populations; we calculated 
these divisions following § 5 of Bell etal. (2004a): 

(t/-V) re st= 1.15-0.3U-0.08(My+20.7). (1) 

Galaxies having rest-frame U—V less than the values provided 
by equation 1 are often referred to as "blue-cloud" galaxies, 
while those with rest-frame U—V greater than the division are 
called "red-sequence" galaxies. The redshift-dependence of 
the blue-cloud/red-sequence galaxy division is thought to be 
largely due to the evolution of the mean age and dust content 
of the blue-cloud population, with a smaller contribution from 
changes in metallicity. In total, we found 2502 blue-cloud 
galaxies and 757 red-sequence galaxies. 

Generally, the blue-cloud and red-sequence populations are 
composed of late-type and early-type galaxies, respectively 
(see, e.g., Bell etal. 2004b, 2005; Mcintosh etal. 2005). 
In support of this point, we have created Figure 2, which 
shows the fraction of E-CDF-S galaxies in our sample hav- 
ing blue-cloud (dashed histogram) and red-sequence (dot- 
ted histogram) colors as a function of galaxy Sersic index 
n. We utilized the Haussler etal. (2007) Sersic indices, 
which were computed using the GEMS zsso images and the 
GALFIT (Peng etal. 2002) two-dimensional light-profile fit- 
ting program. Light-profile studies of large galaxy samples 
have found empirically that a Sersic cutoff of n = 2.5 can 
roughly discriminate between late-type and early-type galax- 
ies (e.g., Blanton etal. 2003; Shen etal. 2003; Hogg etal. 
2004). Galaxies with n £ 2.5 are generally late-type galaxies, 
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FIG. 2. — Distribution of Sersic indices (n) for our initial sample of 
z = 0-1 galaxies in the E-CDF-S having blue-cloud (dashed histogram) and 
red-sequence (dotted histogram) colors. The solid vertical line indicates 
n = 2.5, the empirical cutoff between late-type (n ^ 2.5) and early-type galax- 
ies (n Z 2.5). 

while the majority of galaxies with n ^ 2.5 are early-types 
(vertical line in Fig. 2). Figure 2 shows that there is reason- 
able agreement between late-type and early-type galaxy pop- 
ulations selected using rest-frame optical colors and Sersic in- 
dices. 

We refined further our division between late-type and early- 
type morphologies by visually inspecting the z^so-band im- 
ages of our entire sample of 3259 galaxies to see if there were 
obvious cases where the rest-frame optical colors provided an 
inaccurate morphological classification. For example, late- 
type galaxies that are highly inclined to our line of sight may 
experience significant reddening of the young (and blue) disk 
population and will have red-sequence colors. Also, due to 
variations of the stellar populations in galaxies of a given mor- 
phological class, there will be some scatter in rest-frame U—V 
color near the division of the red-sequence and blue-cloud 
regimes. This will lead to a number of "green" galaxies that 
are misclassified morphologically when the classification is 
based solely on rest-frame color. We stress that such a visual 
inspection is only to correct source classifications that were 
obviously misclassified by simply using color-magnitude di- 
agrams, and our resulting sample should not be considered a 
morphologically-selected sample. 

Based on our visual inspection, we found 76 obvious early- 
type galaxies with blue-cloud colors, and 142 obvious late- 
type galaxies with red-sequence colors; we reclassified these 
sources as late-type and early-type galaxies, respectively. We 
inspected histograms of the fraction of blue-cloud and red- 
sequence galaxies that were reclassified as a function of red- 
shift and find no obvious redshift-related trend. In Fig- 
ure 3, we show zsso-band postage-stamp images of ten ob- 
vious early-type galaxies with blue-cloud colors (top panels), 
as well as ten obvious late-type galaxies with red-sequence 
colors (bottom panels). After reclassifying these objects, we 
were left with 2568 late-type galaxies and 691 early-type 
galaxies with z = 0-1.4. Since we were interested in studying 
the properties of the late-type galaxy populations, we here- 
after refer to our sample of 2568 late-type galaxies as our main 
sample, which we use in subsequent analyses. 

3. PHYSICAL PROPERTIES AND REDSHIFT EVOLUTION OF 
LATE-TYPE GALAXIES 



X-ray Study of Mass-Dependent Star-Formation History of Late-Type Galaxies 
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FIG. 3. — //ST Advanced Camera for Surveys (ACS) Z850 example images often blue-cloud galaxies that have obvious early-type morphologies {upper panel) 
and ten red-sequence galaxies that have obvious late-type morphologies {lower panel) ; these images demonstrate the aid of the HST coverage. Each postage- 
stamp image has a physical size of 30 kpc X 30 kpc, and in each image we indicate the source name {top), the rest-frame U—V color and redshift of the source 
{lower right), and a vertical bar of length 1"5 for scaling reference {lower left); sources harboring an X-ray-detected AGN (see § 4.1) have been annotated with 
the label "AGN" in the lower right. 



The primary goal of this study is to investigate the X-ray 
evolution of normal (i.e., non-AGN) late-type galaxies and 
to determine how this evolution depends upon three intrinsic 
physical properties: optical luminosity, stellar mass, and star- 
formation rate. Below, we describe how we estimated each 
of these physical properties for the galaxies in our main sam- 
ple. We note that the populations of late-type galaxies that 
we are investigating here have been selected via their intrin- 
sic physical properties, which may have changed significantly 
from the observed epoch to the present day. For example, it is 
expected that a significant fraction of the late-type galaxies in 
our main sample will evolve into early-type galaxies via merg- 
ers or passive stellar evolution (e.g., Bell et al. 2007). Ideally, 
we would like to study the evolution of the X-ray properties 



of late-type galaxies while controlling for such changes in the 
physical nature of each galaxy. However, since the details of 
this evolution are highly complex and not well understood for 
a given galaxy, such a task is beyond the scope of this paper. 
We therefore investigate the X-ray evolution of normal late- 
type galaxy populations in relation to their observed intrinsic 
physical properties. 

3.1. Optical Luminosity 

In order to study late-type galaxy samples selected using 
an observable quantity, we made use of the B-band luminos- 
ity (Lg)- In § 2.3, we computed absolute B-band magnitudes 
Mb for galaxies in our main sample using photometrically- 
derived SEDs. The B-band emission from a given late-type 
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galaxy will be significantly influenced by large populations 
of old ( > 100 Myr) stars (measured by the stellar mass) as 
well as the younger and less-numerous massive stars that re- 
side in low-obscuration star-forming regions (measured by the 
star-formation rate). In Figure 4a, we show L B (expressed in 
solar units; where L B q = 5.2 x 10 32 ergs s" 1 ) versus redshift 
for galaxies in our main sample. The redshift-dependent se- 
lection limit of L B for our sample is set by our zsso < 23 crite- 
rion, and at z = 0.5 (z = 1 .4) this limit corresponds to roughly 
L S Ri3x 10 9 Lfl !0 (L B w 4 x 10 10 L B}Q ). For comparison, in 
Figure 4a, we have plotted the values of L B for the Milky 
Way (MW) and local galaxies M101, M82, and the lumi- 
nous infrared galaxy (LIRG) NGC 3256. For the MW, we 
adopted L B = 1.9 x 10 10 L^© as the approximate total disk- 
plus-bulge B-band luminosity (see Table 15.2 of Gilmore et al. 
1990). For the local galaxies, we utilized optical photometry 
from the Third Reference Catalog of Bright Galaxies (RC3; 
de Vaucouleurs etal. 1991) and adopted distances from the 
IRAS Revised Galaxy Sample (Sanders etal. 2003), which 
gives distances of 6.7, 3.6, and 35.4 Mpc and implies corre- 
sponding B-band luminosities of 3.3 x 10 10 , 3.8 x 10 9 , and 
4.8 x 10 10 L BtQ for M101, M82, and NGC 3256, respectively. 

3.2. Stellar Mass 

As discussed in § 1, the X-ray emission from LMXB pop- 
ulations scales with galaxy stellar mass. It is therefore use- 
ful to select late-type galaxies via their stellar masses as a 
means for estimating the LMXB contribution to their X-ray 
emission. To estimate the stellar mass (M*) of each of our 
galaxies, we exploited the tight correlation between rest- 
frame optical color and stellar mass-to-light ratio (e.g., Bell 
& de Jong 2001; Bell etal. 2003; Kauffmann etal. 2003a; 
see also, Borch et al. 2006). For this calculation, we used a 
combination of rest-frame B-V colors and rest-frame TT-band 
luminosities. Zf-band luminosities were computed by fitting 
all available optical/near-IR photometric data to a grid of 276 
synthetic spectra generated by the PEGASE stellar population 
synthesis code (Fioc & Rocca-Volmerange 1997). These tem- 
plates assume a single formation epoch with an exponentially 
decaying star-formation history (time constants r = 1, 2, 4, 
and 7 Gyr) and a Kroupa etal. (1993) IMF; the template grid 
spans ages of 0-15 Gyr. For each galaxy in our main sample, 
we convolved the best-fit template spectrum with the /f-band 
filter function to approximate the rest-frame TT-band luminos- 
ity, L K . This method allowed us to account for significant cur- 
vature in the rest-frame near-IR SEDs of our sources where 
photometric data points did not always overlap with the rest- 
frame #-band. We note that 2546 («99.1%) of our late-type 
galaxies have at least one photometric data point lying at rest- 
frame wavelengths between 1 and 6 /im, which significantly 
constrains the SED in the rest-frame TT-band. We adopted the 
prescription outlined in Appendix 2 of Bell et al. (2003) to 
estimate M* using rest-frame B-V color and TT-band lumi- 
nosity: 

logM„/M Q = logL jC /L^ Q + 0.135(B-y)-0.306. (2) 

The numerical constants in equation 2 were supplied by Ta- 
ble 7 of Bell etal. (2003) and are appropriate for our choice 
of B-V color and L K ; the normalization has been adjusted 
by -0.1 dex to account for our adopted Kroupa (2001) IMF 
(see § 1). We compared our stellar-mass estimates with those 
computed by Borch et al. (2006) for 1758 («68% of our main 
sample) galaxies in the COMBO- 17 survey and find excel- 
lent agreement between methods, with an overall scatter of 




FIG. 4. — Rest-frame B-band luminosity Lg (a), stellar mass (b), and 
star-formation rate SFR (c) versus redshift for sources in our main sample. 
For comparison, in each panel we have indicated the properties of the MW 
and local galaxies M101, M82, and NGC 3256. In Figure 4c, we have in- 
cluded only SFRs for the 888 late-type galaxies in our main sample that have 
24/im counterparts (see § 3.3 for details); the dashed curve indicates roughly 
the SFR detection limit for sources without 24/im counterparts. X-ray- 
detected normal galaxies and AGNs have been indicated with open circles 
and diamonds, respectively. The thick gray rectangles indicate regions where 
galaxy populations were selected for X-ray stacking (see §§4.2 and 5.1). For 
each galaxy. M* and SFR were computed following equations 2 and 3, re- 
spectively. 

w0.2 dex. 

In Figure 4b, we show M + versus redshift for our sample of 
late-type galaxies. We note that our stellar-mass estimates are 
broadly limited by the zs50 < 23 criterion used in our sample 
selection. At z = 0.5 (z = 1.4) our samples are representative 
for late-type galaxies withM* > 1O 9 M (M* > 10 10 M o ). For 
comparison, in Figure 4b we have plotted stellar masses for 
the MW and local galaxies M101, M82, and NGC 3256. We 
adopted a stellar mass of 5 x 10 10 M for the MW (e.g., Ham- 
mer et al. 2007). For the local galaxies, we estimated using 
a similar relation to equation 2, except using B-V color and 
B-band luminosity (see Appendix 2 of Bell et al. 2003 for de- 
tails); for this computation, we utilized B-V colors from RC3 
and values of L B determined in § 3.1. This computation gives 
stellar masses of 1.3 x 10 10 , 5.4 x 10 9 , and 7.8 x 1O 1O M for 
M101, M82, and NGC 3256, respectively. 

3.3. Star-Formation Rates 

Since the X-ray emission from normal late-type galaxies 
is known to be strongly correlated with SFR, it is of partic- 
ular interest to understand how changes in SFR have con- 
tributed to the X-ray evolution of the normal late-type galaxy 
population. To calculate SFRs for the galaxies in our main 
sample, we utilized estimates of both the dust-uncorrected ul- 
traviolet luminosities (Luv) originating from young stars and 
the infrared luminosities (8-1000/im; Lir) from the dust that 
obscures UV light in star-forming regions (see, e.g., Ken- 
nicutt 1998 for a review). The former quantity was com- 
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puted following L uv = 3.3^/^(2800 A), where /„(2800 A) is 
the rest-frame 2800 A monochromatic luminosity (see § 3.2 
of Bell etal. 2005). ^(2800 A) was approximated using 
the template SEDs that were constructed in § 3.2. The lat- 
ter quantity (Lir) was computed using observed-frame 24/im 
flux densities (i.e., rest-frame 24/im/l + z) from observations 
with the MIPS (Rieke etal. 2004) camera onboard Spitzer. 
We note that generally the dust-obscured star-formation ac- 
tivity, probed here using Lir, can be measured using either 
the UV spectral slope or optical nebular recombination lines 
(e.g., Balmer emission lines and O II); however, the data avail- 
able for our z ~ 0-1.4 late-type galaxies lack the rest-frame 
UV and spectroscopic coverage needed to provide such useful 
measurements. 

Over the CDFs, deep 24/im observations were available for 
the GOODS fields (/24^ m , lim ~ 30 /iJy, 6er). 13 For the remain- 
ing area (covering the outer regions of the E-CDF-S), shal- 
lower 24/im observations (/24/jm, lim ~ 120 /iJy) were avail- 
able through SWIRE (Lonsdale etal. 2003). We matched the 
positions of the 2568 late-type galaxies in our main sample 
with those from the available 24/im source catalogs, requir- 
ing that the HST and Spitzer centroids be offset by no more 
than 1"5. We found successful matches for 888 («35%) of 
our late-type galaxies. For these sources, we converted 24/im 
flux densities to Lir following the methods outlined in Pa- 
povich & Bell (2002). Briefly, we utilized the entire grid of 
64 infrared SEDs provided by Dale etal. (2001) to estimate 
the mean conversion factor r)(z), which transforms observed- 
frame 24/im luminosity vl v (24fim/ 1 +z) to Lir as a function 
of redshift [i.e., Lir = 7y(z)zW„(24/im/l +z)]. The mean con- 
version factor spans the tight range of r\ = 7.2-12.2 over the 
entire redshift range z = 0-1.4. We note that different choices 
of infrared SEDs yield similar results at z£ 1, but become 
significantly discrepant at z 1 (see, e.g., Fig. 2 of Papovich 
& Bell 2002). 

Using our estimates of Luv and Lir, we calculated star- 
formation rates for galaxies in our sample using the following 
equation: 

SFR(M Q yr" I ) = 9.8 x 10- n (LiR + L uv ), (3) 

where Lir and Luv are expressed in units of the solar bolo- 
metric luminosity (L Q = 3.9 x 10 33 ergs s" 1 ). Equation 3 
was adopted from § 3.2 of Bell etal. (2005) and was de- 
rived using PEGASE stellar-population models, which as- 
sumed a 100 Myr old population with constant SFR and a 
Kroupa (2001) IMF (see Bell 2003 for further details). We 
note that for the majority of our galaxies Lir > Luv, and for 
sources with z < 0.5 and z > 0.5, the median IR-to-UV ratio 
Lir/Luv = 2.8 and 7.1, respectively. Ourz < 0.5 galaxies have 
median Lir = 1.1 x 10 I0 L q , which is characteristic of objects 
that are found in the nearby D < 10 Mpc universe. By con- 
trast, our z > 0.5 galaxies have median Lir = 9.3 x 1O 1O L0, 
similar to the local LIRG population found at D > 10 Mpc. 
In Figure 4c, we show the distribution of SFRs for the 888 
galaxies in our main sample that had 24/im counterparts. We 
note that for sources in our main sample that were within 
the GOODS regions, where the 24/im observations are most 
sensitive, the infrared detection fraction drops from wl00% 
for galaxies with zsso = 20 ± 0.2 to «20% for galaxies with 
Z830 = 22.8 ± 0.2. This demonstrates that our SFR complete- 
ness is limited primarily by the 24/im sensitivity limit and 

13 Available at http://ssc.spitzer.caltech.edu/legacy/goodshistory.html 



that our sample of zzso < 23 late-type galaxies is highly rep- 
resentative of galaxy populations above the apparent redshift- 
dependent SFR limit shown in Figure 4c (dashed curve). At 
z = 0.5 (z = 1.4), this limit corresponds to SFR « 1 M Q yr -1 
(SFR w 15 Mq yr" 1 ). For comparison, in Figure 4c we 
have plotted the SFRs for the MW and local galaxies M101, 
M82, and NGC 3256. The SFR of the MW was taken to be 
»2 M Q yr" 1 (McKee & Williams 1997), and the SFRs of the 
local galaxies were computed using equation 3. Values of 
/„(2800 A) were approximated from SEDs available through 
the NASA/IPAC Extragalactic Database (NED) 14 andLi R was 
taken from the IRAS Revised Bright Galaxy Sample (Sanders 
etal. 2003). We find SFR = 2.5, 6.2, and 32.2 M yr" 1 for 
M101, M82, and NGC 3256, respectively. 

As a consistency check on our UV-plus-infrared SFR es- 
timates SFR(UV+IR), we calculated radio-derived SFRs us- 
ing 1.4 GHz observations, SFR(1.4 GHz), following equa- 
tion 7 of Schmitt et al. (2006), which we adjusted to be con- 
sistent with our adopted Kroupa (2001) IMF. We matched 
sources in our main sample to 1.4 GHz catalogs, which 
were derived from observations using the VLA in the CDF-N 
(w30 /iJy; Richards etal. 1998) and the ATCA in the 
E-CDF-S («60 /iJy; Afonso etal. 2006; Rovilos etal. 2007). 
Using a matching radius of 1"5, we found a total of 54 
radio sources coincident with our late-type galaxies. We 
found that 48 («89%) of the radio-detected sources had 
24/im counterparts, allowing for reasonable comparison be- 
tween derived SFRs. For these sources we found rea- 
sonable agreement between SFRs derived from UV-plus- 
infrared and radio measurements, with a mean ratio of 
SFR(1.4 GHz)/SFR(UV+IR) = 1.5 ± 1.0. A large number 
of late-type galaxies (840 sources) in our main sample were 
detected in the 24/im observations that were not detected at 
1.4 GHz. All but three of these sources had SFR(1.4 GHz) 
upper limits that were consistent with that expected from es- 
timates of SFR(UV+IR) (see § 4. 1 for further details). 

4. ANALYSIS 

4.1. X-ray-Detected Late-Type Galaxies and AGN 
Identification 

We utilized the multiwavelength observations in the CDFs 
to obtain a census of the active galaxies in our main sample. 
We began by matching the optical positions of our galaxies to 
the X-ray positions of point sources in the CDF catalogs of 
Alexander etal. (2003) 15 for the w2 Ms CDF-N and «1 Ms 
CDF-S and Lehmeretal. (2005b) 16 for the «250ks E-CDF-S. 
For a successful match, we required that the optical and X-ray 
centroids be displaced by no more than 1 .5 times the radius of 
the Chandra positional error circles (80%-90% confidence), 
which are provided in each respective catalog. We note that 
for a small number of the galaxies in our sample at z ^ 0.3, 
moderately luminous off-nuclear X-ray sources (e.g., ultralu- 
minous X-ray sources [ULXs]) that are intrinsically related to 
the galaxies may lie outside of our adopted matching radius; 
we utilized the off-nuclear X-ray source catalog of Lehmer 
etal. (2006) to identify such galaxies and assign X-ray prop- 
erties (see details below). 

14 Available at http://nedwww.ipac.caltech.edu/ 

15 See http://www.astro.psu.edu/users/niel/hdf/hdf-chandra.html for the 
relevant source catalogs and data products for the «2 Ms CDF-N and fslMs 
CDF-S. 

16 See http://www.astro.psu.edu/users/niel/ecdfs/ecdfs-chandra.html for 
the relevant source catalogs and data products for the R3250 ks E-CDF-S. 
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FIG. 5.— Count-rate ratio in the 2-8 keV to 0.5-2 keV band- 
passes (<E>2-8 kev/®0.5-2 kev) versus the logarithm of the 0.5-8 keV flux 
(log/o.5-8 kev) f°r X-ray-detected sources in our main sample. Off-nuclear 
X-ray sources catalogued by Lehmer et al. (2006) have been highlighted with 
open squares. The shaded region represents sources with effective photon 
indices r,.ff £ 1 ; we classified these sources as AGN candidates (see § 4, cri- 
terion 1). For reference, we have plotted lines corresponding to r cff = 0.5, 1, 
and 2 (dashed lines). 



FIG. 6. — S-band magnitude versus log/0.5-8 keV for X-ray-detected 
sources in our main sample. Open circles represent sources that were classi- 
fied as AGN candidates by criterion 1 (see § 4.1 and Fig. 5), and filled circles 
represent all other sources; open squares have the same meaning as in Fig- 
ure 5. The shaded area represents the region where log(/o.5-8 kev///;) > — 1; 
we classified sources in this region as AGN candidates (see § 4.1, criterion 2). 
For reference, the dashed lines represent log(/o.5-8 ^y/ fa) = -2, -1, and 1. 



The Chandra source catalogs were generated using 
wavdetect (Freeman etal. 2002) with false-positive prob- 
ability thresholds of 1 x 10" 7 and 1 x 10" 6 for the Alexan- 
der etal. (2003) and Lehmer etal. (2005b) point-source cat- 
alogs, respectively. However, as demonstrated in § 3.4.2 of 
Alexander etal. (2003) and § 3.3.2 of Lehmer etal. (2005b), 
legitimate lower significance X-ray sources, detected by run- 
ning wavdetect at a false-positive probability threshold of 
1 x 10" 5 , can be isolated by matching with relatively bright 
optical sources; therefore, when matching our late-type galax- 
ies to X-ray-detected sources, we utilized this technique. The 
sky surface density for late-type galaxies in our main sample 
ranges from w23000 deg -2 in the CDF-N to «8800 deg -2 in 
the w250 ks E-CDF-S. The large difference between these 
source densities is primarily due to differences in applied red- 
shift cuts (see § 2.2). Given the fact that the positional un- 
certainties are generally small ( <, 1 "5) for sources within 
6.'0 of the Chandra aim points, as is the case for sources 
in our main sample, the corresponding estimated number of 
spurious matches is small. We estimate that when using 
wavdetect with a false-positive probability threshold of 
1 x 10~ 5 to search for sources in three Chandra bandpasses 
(0.5-2 keV, 2-8 keV, and 0.5-8 keV), we expect sd.8 spuri- 
ous matches. When including the off-nuclear sources from 
Lehmer etal. (2006), we expect an additional w0.5 false 
sources; this brings our total spurious matching estimate to 
«2.3 sources for our main sample. 

Using the matching criteria above, we find that 225 late- 
type galaxies are detected in at least one of the 0.5-2 keV, 
2-8 keV, or 0.5-8 keV bandpasses. Out of these 225 galaxies, 
12 are known off-nuclear X-ray sources from Lehmer etal. 
(2006). Since only one off-nuclear source is detected for each 
host galaxy, we assume that the off-nuclear point-source dom- 
inates the total X-ray emission from each host galaxy. We 
therefore adopted the X-ray properties presented in Table 1 of 
Lehmer et al. (2006) for each off-nuclear host galaxy. 

The unprecedented depths of the CDFs allow for the indi- 
vidual X-ray detection of luminous normal galaxies over the 
entire redshift range of our main sample (z = 0-1 .4); however, 



the majority of the X-ray-detected sources in even the most- 
sensitive regions of the CDFs will be distant (z £ 0.5) AGNs 
(e.g., Bauer et al. 2004), which we want to separate from our 
normal late-type galaxy sample. We identified AGN candi- 
dates using four primary criteria, which utilize (1) X-ray hard- 
ness to identify luminous obscured sources, (2) X-ray-to- 
optical flux ratios to identify additional relatively-unobscured 
AGNs, (3) the X-ray-SFR correlation to identify additional 
lower-luminosity AGNs that are significantly influencing the 
total X-ray emission, and (4) the combination of radio and in- 
frared properties to identify additional Compton-thick AGNs 
and radio-loud AGNs. As a final check on our AGN identifi- 
cations we utilized optical spectroscopic information to iden- 
tify sources with obvious AGN signatures (see criterion 5 be- 
low). In the sections below we provide details of each criteria. 

1. X-ray Hardness. — One unique signature of moderately 
obscured (Nu ~ 10 22 cm" 2 ) AGN activity is a hard X-ray 
spectrum. For normal galaxies, the collective emission from 
X-ray binaries dominates the total 0.5-8 keV power output. 
On average, these sources have observed power-law X-ray 
SEDs with spectral slopes of T w 1.5-1.7 for LMXBs (e.g., 
Church & Balucinska-Church 2001; Irwin etal. 2003) and 
T w 1-2 for HMXBs and ULXs (e.g., Sasaki et al. 2003; Liu 
& Mirabel 2005; Liu etal. 2006); the presence of a signifi- 
cant hot interstellar gas component will steepen the resulting 
X-ray spectral slope (i.e., produce larger T e ff). To identify ob- 
scured AGNs in our sample effectively, we flagged sources 
having effective photon indices of T e ff £ 1 as AGN candi- 
dates. We determined T e ff using the 2-8 keV to 0.5-2 keV 
hardness ratio $2-8 kev/$o.5-2 kev, where $ is the count rate 
for each bandpass. A few sources have only 0.5-8 keV detec- 
tions. Since these sources were not detected in the 0.5-2 keV 
bandpass, our most sensitive bandpass, there must be a sig- 
nificant contribution from the 2-8 keV bandpass such that 
$2-8 kev/^0.5-2 kev ~ 1 (r e ff ~ 0.9). We therefore classified 
these sources as AGN candidates. In Figure 5, we show the 
$2-8 kev / $0.5-2 kev hardness ratio versus the logarithm of the 
0.5-8 keV flux log/o.5-8 kev for the X-ray-detected sources 
in our main sample. The shaded area highlights the region 
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FIG. 7. — Logarithm of the 0.5-8 keV luminosity Lo.5-8keV versus SFR 
for X-ray-detected sources in our sample that have 24/im counterparts. Open 
circles represent sources that were classified as AGN candidates via crite- 
ria 1 and 2 (see § 4.1 and Figs. 5 and 6), and filled circles represent all 
other sources; open squares have the same meaning as in Figure 5. We have 
shown the X-ray-SFR relations calibrated by Bauer etal. (2002b), Ranalli 
etal. (2003), Gilfanov etal. (2004a), and Persic & Rephaeli (2007); each re- 
spective curve has been annotated in the figure. The shaded area above the 
thick dashed line represents the region where Lo.5_g keV is three times larger 
than its value predicted by Persic & Rephaeli (2007); we classified 0.5-8 keV 
detections that lie in this region as AGN candidates (see § 4.1, criterion 3). 



corresponding to T e ff ^ 1 . Sources that had only upper lim- 
its on $2-8 keV / $0.5-2 keV that lie in the shaded region were 
not classified as AGN candidates, while those with lower lim- 
its on $2-8 kev/^o.s-a keV that were in the shaded region were 
flagged as likely AGNs. We note that occasionally an X-ray 
luminous ULX that is too close to the nucleus of its host 
galaxy to be identified as off-nuclear may have T e ff < 1 and 
therefore be classified as an AGN candidate via this criterion; 
however, since only one identified off-nuclear X-ray source 
has r e ff < 1 (open squares in Fig. 5; n.b., upper limits are not 
included), we do not expect such rare sources to significantly 
affect our results. Using criterion 1, we identified a total of 71 
obscured AGN candidates. 

2. X-ray-to-Optical Flux Ratio. — Detailed analyses of the 
X-ray spectra of luminous AGNs in the « 1 Ms CDF-S show 
that the intrinsic AGN power-law photon index is relatively 
steep, (r) = 1.75 ±0.02 (e.g., Tozzi etal. 2006). Therefore, 
luminous AGNs having column densities of Nh ^ 10 22 cm" 2 
will often have effective photon indices of T e ff > 1 and would 
not have been classified as potential AGNs by criterion 1. 
In order to identify luminous AGNs with Nh ^ 10 22 cm" 2 , 
we utilized the X-ray-to-optical flux ratio (/0.5-8 keV / fid as 
a discriminator of AGN activity (e.g., Maccacaro etal. 1988; 
Hornschemeier etal. 2000; Bauer etal. 2004). We identified 
sources with log(/b.5-8 keV //«) > -1 (see criterion 3 below for 
justification) as unobscured AGN candidates. In Figure 6, we 
show the /?-band magnitude versus log /0.5-8 keV for sources 
in our main sample; the shaded area shows the region where 
l°g(/o.5-8 kev//fi) > -1- Sources that were classified as AGN 
candidates via criterion 1 are denoted with open symbols. In 
total, 43 X-ray-detected sources satisfied criterion 2, and 16 
of these sources were uniquely identified using this criterion 
(i.e., not identified by criterion 1). 

3. X-ray-to-SFR Correlation. — Taken together, criteria 1 
and 2 provide an effective means for identifying AGNs that 
are affected by large absorption column densities (criterion 1) 



and those that are notably X-ray overluminous for a given 
optical luminosity (criterion 2). However, these criteria will 
still miss moderately luminous unobscured AGNs that have 
l°g(/o.5-8 kev//«) < -1 (see, e.g., Peterson etal. 2006). Al- 
though an accurate classification for all such sources is cur- 
rently not possible, the situation can be mitigated using the 
available multiwavelength data. We therefore exploited the 
correlation between X-ray luminosity L x and SFR (see § 1) to 
identify additional AGN candidates in our main sample that 
have significant X-ray excesses over what is expected based 
on the Lx-SFR correlation. In order to calculate the rest-frame 
luminosity (L El - El ; where E\ and £2 are the photon-energy 
lower and upper bounds, respectively) of a source having a 
power-law SED, we used the following equation: 

L El . E2 =4nd 2 L f Ei - E2 (l+z) r - 2 , (4) 

where f Ei - El is the observed-frame emission in the E\-E2 
bandpass and d L is the luminosity distance. Using equa- 
tion 4 and an adopted photon index of T = 2, we calcu- 
lated 0.5-8 keV luminosities for X-ray-detected sources in 
our sample. In Figure 7, we show the logarithm of the 
0.5-8 keV luminosity logLo.5-8 keV versus SFR (computed 
following equation 3) for galaxies in our main sample that 
had 24^m counterparts (see § 3.3). Several estimates of the 
L X -SFR correlation have been shown for reference (Bauer 
etal. 2002a; Ranalli etal. 2003; Gilfanov etal. 2004a; Per- 
sic & Rephaeli 2007, hereafter PR07); these correlations have 
been corrected for differences in X-ray bandpass and SED as 
well as adopted IMFs. Hereafter, we adopt the L X -SFR corre- 
lation from PR07 for comparisons; however, the use of other 
Lx-SFR correlations would yield similar results and conclu- 
sions. Open squares show the locations of the galaxies hosting 
off-nuclear sources, which appear to be preferentially located 
near the L X -SFR correlation. Open circles indicate sources 
that were identified as AGN candidates via criteria 1 and 2. 
Generally, these AGNs have Lo.5-8 kev/SFR ^ 3 times that 
predicted by the PR07 Lx-SFR correlation (thick dashed line), 
which corresponds to a factor of ^ 2.5 times the RMS scat- 
ter of the PR07 L X -SFR correlation. We therefore classified 
all sources in this regime (shaded region) as AGN candidates. 
We note that sources having only 0.5-8 keV upper limits were 
not classified as AGNs. However, sources that were detected 
in the 0.5-8 keV bandpass that had only upper limits on the 
SFR were classified as AGN candidates if they were within 
the shaded region of Figure 7. Using this criterion, we identi- 
fied 101 potential AGN candidates, of which 33 were unique 
to criterion 3. For the 124 X-ray-detected sources that were 
not classified as AGNs, we found that 19 sources («15%) had 
X-ray spectral properties indicative of AGNs (criterion 1) and 
only one source (w0.8%) had an X-ray-to-optical flux ratio 
elevated in comparison to normal galaxies (criterion 2). 

4. Infrared and Radio Properties. — In addition to the 
X-ray-detected sources, we also expect there to be additional 
AGNs that lie below the X-ray detection threshold that are not 
identified here. In § 3.3, we noted using 48 sources that the 
SFRs derived from the UV-plus-infrared emission are gener- 
ally consistent with those derived from radio ( 1 .4 GHz) obser- 
vations; however, there are a few sources that are significantly 
scattered (i.e., by more than a factor of five times the intrin- 
sic scatter) outside of the SFR(UV+IR)-SFR(1.4 GHz) cor- 
relation. These sources have either (1) excess infrared emis- 
sion due to reprocessed dust emission from a highly-obscured 
(Nh ~ 10 24 cm" 2 ) moderately luminous (intrinsic Lx ~ 10 42 - 
10 43 ergs s" 1 ) Compton-thick AGN (e.g., Daddi et al. 2007a,b) 
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FIG. 8. — Logarithm of the 0.5-8 keV luminosity Lo.5-8 keV versus redshift 
for the 225 X-ray-detected sources in our main sample. Open circles indicate 
all of the AGN candidates that were isolated using the four criteria outlined 
in § 4.1, and filled circles represent potential normal galaxies; open squares 
have the same meaning as in Figure 5. The estimated 0.5-8 keV detection 
limits for the «2 Ms CDF-N, the »1 Ms CDF-S, and the «250 ks E-CDF-S 
have been indicated with dotted, dashed, and dot-dashed curves, respectively. 

or (2) an excess of radio emission due to prominent radio- 
emitting jets (e.g., Snellen & Best 2001). The reasonable 
agreement between SFR(UV+IR) and SFR(1.4 GHz) suggests 
that such AGNs are not prevalent in our late-type galaxy sam- 
ple. Ideally we would also like to explore how UV-corrected 
SFRs compare to the UV-plus-infrared SFR to search for po- 
tential Compton-thick AGNs following the approach used by 
Daddi et al. (2007 a,b) for BzK galaxies at z £ 2; however, 
since our galaxies generally lie at much lower redshifts that 
the Daddi et al. (2007a,b) sources, we are unable to constrain 
well UV spectral slopes with our available photometry, which 
would be required to make dust-corrections to the observed 
UV SFRs (see § 3.3 for discussion). For our sample of 24/im- 
detected galaxies, we identified three radio-excess galaxies 
with SFR(1.4 GHz)/SFR(UV+IR) > 5 and three IR-luminous 
galaxies with SFR(1.4 GHz)/SFR(UV+IR) < 1/5. Using 
these criteria we identified three new AGN candidates, two ra- 
dio luminous (J33 146.6-275735 and J33244.3-275141) and 
one IR luminous (J3 3240.0-2742 14); one of these sources, 
the radio luminous source J33244. 3-275 141, was detected in 
the 0.5-2 keV and 0.5-8 keV bandpasses. 

As an additional test, we utilized IRAC photometry (see 
§ 2.3) to search for infrared power-law sources having near- 
IR spectral properties characteristic of luminous AGNs (e.g., 
Alonso-Herrero etal. 2006; Donley etal. 2007). When 
searching for power-law sources, we adopted the criteria dis- 
cussed in Donley et al. (2007). We found that no sources in 
our sample satisfied these criteria, which is consistent with the 
finding that most IRAC power-law sources reside at z ^ 1 ■ 

5. Optical Spectroscopy. — As a final check on our AGN 
classifications, we searched the optical spectroscopic catalogs 
available for CDF sources (e.g., Barger etal. 2003; Le Fevre 
etal. 2004; Szokoly etal. 2004; Wirth etal. 2004; Vanzella 
etal. 2005, 2006) to isolate additional luminous AGNs in our 
sample. In total, we found 15 galaxies in our main sam- 
ple that were classified as AGNs via optical spectroscopy 
and all of these sources had been identified as AGN candi- 
dates by the previous criteria. We note that the majority of 
the X-ray-detected AGNs have moderate luminosities (intrin- 
sic L x ~ 10 41 — 10 43 ergs s" 1 ) and therefore often have high- 



excitation AGN emission lines that are too faint with respect 
to stellar emission to be identified via optical spectroscopy 
(e.g., Moran etal. 2002). Therefore, it is not surprising that 
we do not find any additional AGNs using this criterion. 

To summarize, in the X-ray band we have detected a total 
of 225 (~9%) late-type galaxies out of the 2568 sources in 
our main sample. Using the criteria presented above, we clas- 
sified 121 X-ray-detected sources as AGN candidates (with 
an additional two X-ray-undetected AGN candidates via cri- 
terion 4). The remaining 104 X-ray-detected sources that we 
do not classify as AGN candidates are considered to be nor- 
mal late-type galaxies, and we include these galaxies in sub- 
sequent X-ray stacking analyses (see details in § 5 below). 
Thus we use 2447 late-type galaxies in our stacking analyses. 
In Figure 8, we show log Lo.5-8 keV versus redshift for X-ray- 
detected sources in our main sample, and in Table 1 (available 
electronically) we summarize their properties. AGN candi- 
dates are denoted with open circles and normal galaxies are 
plotted with filled circles. We note that the above criteria 
are not completely sufficient to classify all X-ray-detected 
sources that are truly AGNs as AGN candidates (see, e.g., 
Peterson etal. 2006). Such a misclassification is possible 
for low-luminosity AGNs that are only detected in the more 
sensitive 0.5-2 keV bandpass, which have 2-8 keV emission 
too weak for an accurate classification. In § 5.2, we use the 
2-8 ke V AGN fraction as a function of X-ray luminosity to ar- 
gue quantitatively that we do not expect misclassified AGNs 
(detected only in the 0.5-2 keV bandpass) and low-luminosity 
AGNs below the X-ray detection limit to have a serious im- 
pact on our results. 

4.2. X-ray Stacking Analyses of Normal Late-Type Galaxy 
Populations 

The majority of the normal late-type galaxies that make up 
our main sample were not detected individually in the X-ray 
bandpass. In order to study the mean X-ray properties of 
these sources, we implemented stacking analyses of galaxy 
populations selected by their physical properties (see § 3). 
We divided our main sample into subsamples (to be used for 
stacking) of normal late-type galaxies selected by both physi- 
cal properties (i.e., fi-band luminosity, stellar mass, and star- 
formation rate) and redshifts. In Figure 4, we have highlighted 
the divisions of our sample with thick gray rectangles, and for 
normal late-type galaxies in each region, we used X-ray stack- 
ing analyses to constrain average properties. 

For each of the subsamples defined above, we performed 
X-ray stacking in each of the three standard bandpasses (i.e., 
SB, HB, and FB; see § 1). We expect these bandpasses 
to sample effectively power-law X-ray emission originating 
from X-ray binaries (i.e., HMXBs and LMXBs) with a minor 
contribution from hot interstellar gas in the SB for late-type 
galaxies at z£ 0.5. In our analyses, we used data products 
presented in Alexander et al. (2003) for the «2 Ms CDF-N 
and w 1 Ms CDF-S and Lehmer et al. (2005b) for the «250 ks 
E-CDF-S (see footnotes 15 and 16). Our stacking procedure 
itself was similar to that discussed in § 3 of Steffen etal. 
(2007). This procedure differs from past stacking analyses 
(e.g., Lehmer et al. 2007) in how the local X-ray background 
of each stacked sample is determined, and produces results 
that are in good agreement with the method discussed in § 2.2 
of Lehmer etal. (2007). For completeness, we have outlined 
this procedure below. 

Using a circular aperture with radius R dp = l."5, we ex- 
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Source Name zsso L B M* SFR £o.5-8 keV fo.ss keV 

(J2000.0) (mag) z (logL B , ) (logM Q ) (Mq yr 1 ) Survey (ks) (log ergs cm" 2 s" 1 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



J033132.81-280115.9 17.19 0.15P 10.39 10.49 <5.89 E-CDF-S 03 217 -15.06 

J033132.84-280107.5 22.17 0.15? 8.85 8.44 0.23 E-CDF-S 03 217 -15.45 

J033139.20-280222.3 17.06 0.25P 10.94 11.07 1.40 E-CDF-S 03 206 < -15.15 

J033139.97-274157.0 20.14 0.47 s 10.14 10.43 <19.93 E-CDF-S 02 226 -13.46 

J033141.79-275635.2 19.58 0.5F 10.84 10.54 <13.17 E-CDF-S 03 239 -15.02 



NOTE. — Col.(l): Chandra source name. Col. (2): ACS egso-band magnitude. Col. (3): Redshift estimate. Superscripts "s" and "p" indicate spectroscopic and photometric redshifts, 
respectively (see § 2 for details). Col. (4): Logarithm of the rest-frame B-band luminosity in units of Lb.q ■ Col. (5): Logarithm of the stellar mass in units of Mq . Col. (6): Star-formation 
rate in units of Mq yr -1 . Col.(7): Survey field in which each source was identified. For E-CDF-S identifications, the associated field number (i.e., 01-04) indicates the Chandra pointing 
within which the source was detected (see Lehmer et al. 2005b for details). Col. (8): Effective 0.5-8 keV exposure time (in units of ks). Col.(9)-(l 1): Logarithm of the 0.5-8 keV, 
0.5-2 keV, and 2-8 keV flux in units of ergs cnT 2 s~'. Col.(12)-(14): Logarithm of the 0.5-8 keV, 0.5-2 keV, and 2-8 keV rest-frame luminosity in units of ergs s" 1 . Col.(15): 
Effective photon index (r cff ). Col.(16): Logarithm of the 0.5-8.0 keV to R-band flux ratio. Col.(17): AGN candidate (Y/N)? Col.(18): AGN selection criteria used (i.e., 1-5; see § 4.1 
for details). All 225 entries of Table 1 are available in the electronic edition. A portion is shown here for guidance regarding its form and content. 



Table 2. Stacked Late-Type Normal Galaxies: Mean Properties 



Net Source Counts (S-B) Signal-to-Noise (S/N) 



Selection Type z msm /V ga i N iet 0.5-8 keV 0.5-2 keV 2-8 keV 0.5-8 keV 0.5-2 keV 2-8 keV 
(1) (2)_ " 

logLg/Ls © = 9.5-10.0 0.23 ± 0.06 

logL B /.,•. = 9.5-10.0 0.46 ± 0.06 

\ogL B /L B q = 10.0-10.5 0.24 ±0.02 

logL B /L B]Q = 10.0-10.5 0.33±0.03 

UwJn /.a]. = 10.0-10.5 0.41 ±0.03 

logL B /L BQ = 10.0-10.5 11.511 : 0.02 

logL B /L B Q = 10.0—10.5 11.57 : 11.02 

\ogL B /L B Q = 10.0—10.5 0.67±0.02 

logL B /L BQ = 10.0-10.5 0.75 ±0.02 

logZ.fi/LB0 — 10.0—10.5 0.84±0.02 



(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 
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3.0 
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1 
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6.8 


0.9 


28 


2 


66.9 


46.5 


20.2 


5.2 


5.3 


2.1 


33 


3 


82.9 


68.5 


14.9 


5.9 


6.7 


1.5 


67 


7 


204.2 


144.7 


59.8 


9.3 


9.7 


3.8 


158 


6 


347.5 


240.0 


109.1 


10.7 


11.3 


4.4 


100 


3 


128.0 


109.0 


18.7 


6.0 


7.5 


1.2 


71 


2 


154.1 


93.5 


60.6 


6.3 


6.1 


3.2 


62 


2 


156.3 


107.4 


48.8 


6.4 


6.9 


2.6 


57 


3 


105.9 


64.8 


41.2 


4.4 


4.5 


2.2 



NOTE. — Col.(l): Physical property (i.e., L B , M„, and SFR) used to select the stacked sample. Col. (2): Mean redshift (z me an) and 1 a standard deviation. Col. (3): Number of 
galaxies stacked (iVgai). Col. (4): Number of X-ray-detected normal galaxies stacked. Col.(5)-(7): Net source counts (S-B) for the 0.5-8 keV, 0.5-2 keV, and 2-8 keV bandpasses. 
Col.(8)-(10): Signal-to-noise ratio (S/N) for the 0.5-8 keV, 0.5-2 keV, and 2-8 keV bandpasses. Col.(l 1)— (13): Logarithm of the mean 0.5-8 keV, 0.5-2 keV, and 2-8 keV flux 
in units of ergs cirT 2 s' 1 . Col.(14)-(16): Fraction of the mean 0.5-8 keV, 0.5-2 keV, and 2-8 keV flux originating from the X-ray-undetected galaxies. Col.(17)-(19): Logarithm 
of the mean 0.5-8 keV, 0.5-2 keV, and 2-8 keV rest-frame luminosity. Col.(20): Mean effective photon index (IV). Col.(21): Logarithm of the 0.5-8.0 keV to S-band flux ratio. 
Col.(22)-(24): Mean values of L B , M„ and SFR for samples selected by L B , M„ and SFR, respectively. Col.(25)-(27): Logarithm of the ratios L X /L B (ergs s' 1 Lj' Q ), L x /M„ (ergs s" 1 

M^J), and Lx/SFR [ergs s' 1 (M e yr 1 ) -1 ] for samples selected by L B , M,, and SFR, respectively. Col.(28)-(30): Estimated fraction of the mean 0.5-8 keV, 0.5-2 keV, and 2-8 keV 
stacked emission originating from undetected AGNs. All 44 entries of Table 2 are available in the electronic edition. A portion is shown here for guidance regarding its form and 
content. 



tracted Chandra source-plus-background counts 5, and expo- 
sure times Tj (in units of cm 2 s) for each galaxy using images 
and exposure maps, respectively. For a given source, we used 
only Chandra pointings with aimpoints (see footnote 12) that 
were offset from the source position by less than 6f0; here- 
after, we refer to this maximum offset as the inclusion radius, 
^inci- 17 Since the Chandra PSF increases in size with off- 
axis angle and degrades the sensitivity for sources that are 
far off-axis, our choices of R. dp and R[ nc i have been chosen 
to give the maximal stacked signal with the majority of the 
PSFs being sampled by our stacking aperture (see § 2.2 and 
Fig. 3 of Lehmer et al. 2007 for justification). We estimate, 
based on stacked images of z ~ 0.1-0.3 optically luminous 
(Lb > 2 x 10 9 Lb,q) late-type galaxies, that £ 10% of the nor- 
mal galaxy X-ray emission originates outside of our 1 "5 ra- 
dius aperture (2.8-6.7 kpc at z = 0.1-0.3) for all stacked sam- 
ples. For galaxies that were within 6.'0 of more than one of the 
CDF aimpoints, we added source counts and exposure times 

17 Note that S ind has the same meaning as it did in Lehmer et al. (2007). 



from all appropriate images and exposure maps, respectively; 
however, there were very few sources in our main sample that 
met this criterion. 

Using background maps (see § 4.2 of Alexander et al. 2003 
for the w2 Ms CDF-N and wl Ms CDF-S and § 4 of Lehmer 
etal. 2005b for the «250 ks E-CDF-S) and exposure maps, 
we measured local backgrounds B,.i 0C ai and exposure times 
Ti.iocai within a 30 pixel x 30 pixel («15" x 15") square, cen- 
tered on each source with the 1 "5 radius circle masked out. 
Here again, if a source was within 6'.0 of more than one of the 
CDF aimpoints, we summed the local backgrounds and expo- 
sure times. We estimated the expected number of background 
counts in each circular aperture B, by scaling the background 
counts within the square by the relative exposure times of the 
circular aperture and the square (i.e., 5, = B^\ oc . d i x 7}/7],i OC ai)- 
This approach is similar to scaling the background counts in 
the square by the relative areas of the circular aperture and 
the square; however, by using exposure times, we are able to 
account more accurately for spatial variations in pixel sensi- 
tivity due to chip gaps, bad pixels, and vignetting. Further- 
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more, comparisons between this method and the Monte Carlo 
method used in Lehmer et al. (2007) to compute B, give ex- 
cellent agreement and are most convergent for large numbers 
of Monte Carlo trials. 

When stacking galaxy populations, we excluded sources 
that were (1) classified as AGN candidates (via the crite- 
ria outlined in § 4.1), (2) within 10"of an unrelated source 
detected in the X-ray catalogs, (3) within the extent of ex- 
tended X-ray sources (see Bauer etal. 2002b, § 3.4 of Gi- 
acconi etal. 2002, and § 6 of Lehmer etal. 2005b), and (4) 
located within 3"of another late-type galaxy in our main sam- 
ple. We note that we include X-ray-detected normal galax- 
ies when stacking our samples, since we are interested in the 
average properties of the normal late-type galaxy population. 
We tested the effects of including such sources by stacking 
samples both with and without X-ray-detected sources in- 
cluded, and find similar results for both cases. Typically, the 
X-ray-detected sources contribute w 10-30% of the counts in 
the total stacked signal. For each of the subsamples of normal 
late-type galaxies outlined in Figure 4 (gray rectangles), we 
determined stacked source-plus-background (S = ^jSf) and 
background counts (B = J2i Bd to determine net counts (S-B). 
For each stacked sample, we required that the signal-to-noise 
ratio [S/N = (S-B)/VB; where \S-B\ <C B and B > 20] be 
greater than or equal to 3 (i.e., ^ 99.9% confidence) for a de- 
tection. For stacked samples without significant detections, 
3cr upper limits were placed on the source counts. 

We converted the net counts obtained from each stacked 
sample to absorption-corrected fluxes and rest-frame lumi- 
nosities using a power-law SED with T = 2. Due to the fact 
that our l."5 radius stacking aperture encircles only a frac- 
tion of the PSF 18 for sources at relatively large off-axis an- 
gle, we calculated aperture corrections £, for each stacked 
source i. Since we are calculating average X-ray counts from 
the summed emission of many sources of differing back- 
grounds and exposure times, we used a single, representative 
exposure-weighted aperture correction, £. This factor, which 
was determined for each stacked sample, was calculated as 
follows: 

(5) 

where T = V. 7}. The average aperture corrections (£) for 
sources in our main sample were wl.6, 1.8, and 1.7 for the 
SB, HB, and FB, respectively. Using our adopted SED, we es- 
timated observed mean X-ray fluxes using the following equa- 
tion: 

fEi-E 2 =A Ei -e 2 £ ( — 7— ) , (6) 



where Ae x -e 2 is a bandpass-dependent count-rate-to-flux con- 
version factor that incorporates both the X-ray SED informa- 
tion as well as Galactic extinction using the column densities 
listed in § 1 . These mean X-ray fluxes were then converted 
to rest-frame luminosities following equation 4, assuming a 
photon index of T = 2. 

5. RESULTS 
5.1. Stacking Results 

18 At off-axis angles 9 f» 3', our 1"5 radius circular aperture contains an 
encircled-energy fraction of sdl00%, 80%, and 100% for the SB, HB, and 
FB, respectively; however at 9 & 6', this fraction decreases to sd30%, 25%, 
and 25%, respecively. 
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FIG. 9. — Effective photon index (r c ff) versus the logarithm of the X-ray- 
to-optical flux ratio (log/o.5-8 kcv/Zs) f° r 44 stacked samples selected via 
observed properties: Lg (filled circles), M k (filled squares), and SFR (filled 
triangles). The characteristic mean error bar for each quantity is given in the 
lower right-hand corner. The median logarithm of the X-ray-to-optical flux 
ratio is indicated with a vertical dotted line (log/0 5-8 kev//i? = -2.53). The 
median effective photon index for the samples that were detected in both the 
SB and HB is indicated with a horizontal dotted line (r™= dian = 1.47). We 
note that all of the stacked samples with both SB and HB detections have 
T e ff > 1, consistent with normal galaxies. For the several stacked samples 
that have only SB detections, we have indicated lower-limits on r cff . The 
shaded regions and corresponding boundaries (dashed lines) represent areas 
where X-ray-detected sources were classified as AGN candidates (for details, 
see discussion of criteria 1 and 2 in § 4.1); note that additional criteria were 
used to identify potential AGNs when generating our samples of normal late- 
type galaxies (see criteria 3-5 in § 4.1). 



Using the stacking analysis methods discussed in § 4.2, we 
stacked the late-type galaxy samples presented in Figure 4 
(stacked samples are denoted with thick gray rectangles). 
These samples were selected using L B , M+, and SFR, which 
include 14, 17, and 13 stacked samples (44 total), respectively. 
In Table 2 (available electronically), we tabulate our X-ray 
stacking results. We found significant (i.e., S/N >, 3) X-ray 
detections in the 0.5-2 keV and 0.5-8 keV bandpasses for all 
stacked samples. In the 2-8 keV bandpass, 15 out of the 44 
stacked samples were detected, and these samples generally 
constitute the most optically luminous and massive galaxies, 
as well as those galaxies that are most actively forming stars. 

In Figure 9, we show the effective photon index (r e ff) ver- 
sus the logarithm of the X-ray-to-optical mean flux ratio 
(log/0.5-8 keV / /r) for our 44 stacked samples (filled symbols) 
that were selected via their observed properties. Effective 
photon indices were estimated using HB-to-SB count-rate ra- 
tios (i.e., $2-8 kev/^o.5-2 kev)- All 44 stacked samples have 
X-ray-to-optical flux ratios and X-ray spectral slopes consis- 
tent with normal galaxies (unshaded region in Fig. 9), sug- 
gesting that these samples are not heavily contaminated by an 
underlying population of AGNs. 19 We find a median loga- 
rithm of the X-ray-to-optical flux ratio of log/b.5-8 kev //« = 
-2.53 (vertical dotted line), and for samples that were de- 
tected in both the HB and SB, the median effective photon 

19 We note that two of our samples have spectral slopes that are near our 
chosen division between normal galaxies and AGNs (r c ff ~ 1). These sam- 
ples have logLg = 10-10.5 and logM* = 9.5-10.1 and lie at z K, 0.66, a red- 
shift where large-scale sheets of galaxies and AGNs have been isolated pre- 
viously in the CDF-S (e.g., Gilli etal. 2003).. These galaxies have spectral 
slopes that are either statistically scattered towards r cff = 1 (see characteris- 
tic error bar in Fig. 9) or may have some non-negligible contribution from 
heavily-obscured AGNs. 
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FIG. 10. — Logarithm of (a) the X-ray-to-B-band mean luminosity ra- 
tio (\ogLx/L B ), (b) the X-ray-to-stellar-mass mean ratio (logLx/Mt), and 
(c) the X-ray-to-star-formation-rate mean ratio (log Lx /SFR) versus redshift 
{filled symbols and curves) for stacked normal late-type galaxy samples se- 
lected by L B (Fig. 4a), M* (Fig. 4b), and SFR (Fig. 4c), respectively. For 
comparison, in each panel we have indicated the properties of the MW and 
local galaxies Ml 01, M82, and NGC 3256. Quoted X-ray luminosities corre- 
spond to the 0.5-8 keV bandpass and were calculated following the methods 
described in § 4.2, assuming a power-law SED with photon index of T = 2. 
Symbols and curves correspond to unique ranges of L B , M t , and SFR, which 
are annotated in each respective figure. For reference, in Figures 10a and 10b 
we have plotted the corresponding values of \og Lx/L B and log Lx/M*, re- 
spectively, for normal late-type galaxies in the local universe (open symbols) 
using the S01 sample. In Figure Wb we show the expected LMXB contribu- 
tion based on Gilfanov etal. (2004b; dashed line). Finally, in Figure 10c we 
show the local Lx-SFR relation and its dispersion (dashed line with shading) 
derived by PR07 and corrected for our choice of IMF and T = 2. We note 
that roughly all of our data points in Figure 10c lie «0. 1-0.2 dex below the 
PR07 relation, which is likely due to systematic differences in how Lx and 
SFR were determined between studies. 



index is r^ dlan = 1.47 {horizontal dotted line). These values 
are characteristic of galaxies dominated by X-ray binary pop- 
ulations. 

In Figure 10, we show the logarithm of the ratio of the 
0.5-8 keV luminosity (hereafter, Lx) to each physical prop- 
erty (i.e., logLx/L B , logLx/AL, and logL x /SFR) versus red- 
shift for our samples. Each relevant quantity for the MW 
and local galaxies M101, M82, and NGC 3256 have been 
shown for comparison. We adopt 6.6 x 10 39 ergs s" 1 as 
the approximate X-ray luminosity of the MW (Grimm et al. 
2002). For the local galaxies, we utilized the X-ray lumi- 
nosities from Shapley etal. (2001; hereafter S01) for M101 
(L x = 1.6 x 10 40 ergs s" 1 ) and M82 (L x = 4.6 x 10 40 ergs s" 1 ) 
and Lira et al. (2002) for NGC 3256 (L x = 4.0 x 10 41 ergs s" 1 ). 
We have corrected the X-ray luminosities to be consistent 
with our use of the 0.5-8 keV bandpass and our choice of 

r = 2. 

For the purpose of comparing our results to those for nearby 
late-type galaxy populations, we made use of the S01 sam- 
ple of 183 normal local (D £ 100 Mpc) spiral and irregular 
galaxies. These galaxies were observed in the 0.2-4 keV band 
using the Einstein IPC and HRI (Fabbiano etal. 1992) and 



AGNs, with luminous X-ray emission or spectral signatures 
indicative of AGN activity, have been excised from the sam- 
ple. To avoid the inclusion of early-type SO galaxies, we chose 
to utilize 139 normal late-type galaxies from the S01 sample 
with morphological types T > 2 (Sa-Irr Hubble types). 

The S01 sample covers ranges of Lb and AL that are well- 
matched to our main sample and are representative of late- 
type galaxies in the local universe. In order to compute mean 
X-ray luminosities for samples that were directly comparable 
with our results, we divided the S01 sample into the same in- 
tervals of Lb and that were used for our main sample (see 
Figs. 4a and 4b). Since several of the S01 galaxies had only 
X-ray upper limits available, we computed mean X-ray lumi- 
nosities and errors using the Kaplan-Meier estimator available 
through the Astronomy SURVival Analysis software pack- 
age (ASURV Rev. 1.2; Isobe & Feigelson 1990; LaValley 
et al. 1992); the Kaplan-Meier estimator handles censored data 
sets appropriately. When computing these mean X-ray lumi- 
nosities, we filtered the S01 samples appropriately into dis- 
tance intervals to avoid the Malmquist bias. In Figures 10a 
and 10b, we show the corresponding values of log Lx/Lb and 
log Lx/M+, respectively, for the S01 sample with open sym- 
bols. By contrast, the SFRs of the local sample are generally 
too low ( £ 1-10 Mq yr" 1 ) to provide a meaningful compar- 
ison with our distant 24,um-detected sources. This is due to 
the strong positive evolution of the star-formation rate density 
with redshift (see § 1), which makes SFRs that are common 
for galaxies in our sample ( ^ 1-10 M Q yr" 1 ) comparatively 
rare at z = 0. 

From Figures 10a and 10b, it is apparent that there is signif- 
icant positive redshift evolution in log Lx/Lb and log Lx/M* 
over the redshift range of z ~ 0-1.4. For each of the six to- 
tal selection ranges of Lb and M*, the redshift progression 
of X-ray luminosities is inconsistent with a constant at the 
>99.9% confidence level. For the most optically luminous 
(L B =3-20 x 10 10 L B)Q ) and massive (Af* = 1-20 x 10 10 M Q ) 
late-type galaxies at z = 1.4, Lx/Lb and Lx/M* are mea- 
sured to be larger than the local values of S01 by factors of 
3.1 ±0.7 and 9.6 ±3.1, respectively. Such values are consis- 
tent with the ss(l +z) L5 " 3 evolution of L x found for the nor- 
mal late-type galaxy population, which has been constrained 
using largely the most optically luminous and massive galax- 
ies (Georgakakis etal. 2006; Ptak etal. 2007; Tzanavaris & 
Georgantopoulos 2008). 

The above results confirm the increase in L X /L B with red- 
shift found by H02. In past studies, Lx/Lb has been used 
as a proxy for star-formation activity (e.g., Ptak etal. 2001; 
H02; Lehmer et al. 2005a), despite the fact that Lb is likely 
to be somewhat sensitive to SFR (see discussion in § 3.1). 
As discussed in § 1, L x has been shown to be strongly cor- 
related with galaxy SFR. Correlation studies of spiral galax- 
ies in the local universe have also found strong correlations 
between Lx and Lb such that Lx oc L B 5 (S01; Fabbiano & 
Shapley 2002). The nonlinear relationship between X-ray 
and B-band emission is thought to be due to the increase in 
dust obscuration with star-formation activity, which attenu- 
ates light from the B-band more effectively than it does in 
the X-ray band. Using our sample of 47 normal late-type 
galaxies with both X-ray and 24[im detections, we found that 
Lx/Lb was positively correlated with the UV dust-extinction 
measure (L IR +L UV )/L UV (Kendall's r = 0.44; > 99.99% con- 
fidence level), thus providing support for this hypothesis. Fur- 
thermore, it has been shown that Lx/Lb is correlated with 
^60/im/Lioo^m, suggesting that more intense X-ray emission 
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Table 3. Parametric Fitting Results For Stacked Samples 



/phys 


V 


x 2 


A 


B 


C 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


L B 


11 


25.37 


31.34±1.23 


0.85 ±0.13 


2.91 ±0.61 


M* 


14 


25.73 


33.97 ±0.87 


0.58 ±0.07 


4.48 ±0.37 


SFR 


10 


9.72 


39.91 ±0.06 


0.89±0.11 


0.32 ±0.88 



NOTE. — This table contains basic fitting parameters for \ 2 nts to our X-ray stacking 
results. For each sample, selected by physical property /phys, we performed parametric 
fits for the mean 0.5-8 keV luminosity Lx following logLx = A+Slog/ p h ys +Clog(l+2). 
Col.(l): Physical parameter / p h ys used in fitting our stacking results. Col. (2): Number 
of degrees of freedom used in each fit. Col. (3): Minimum \ 2 value for each fit. Col. (4)- 
(6): Best-fit values of A, B, and C with errors (90% confidence). For futher details, see 
§5.1. 



is associated with hotter IR colors, which is indicative of in- 
tense obscured star-formation activity. Therefore, the increase 
in L X /L B with redshift (see Fig. 10a) observed for our galax- 
ies is likely due to an increase in their star-formation activity. 

Upon comparing Lb with M* and SFR for our galaxies, 
we find that Lb is well correlated with both (Kendall's 
t = 0.69) and SFR (r = 0.58); however, the scatter in the L B - 
M+ relation is significantly smaller than it is for the L#-SFR 
relation (w0.2 dex versus w0.4 dex, respectively) implying 
that Lb traces more effectively rather than SFR. It is there- 
fore not surprising that we see similar redshift evolution of 
L x /M+ and L X /L B . 

In Figure 10b, we show the estimated LMXB contribution 
to \ogL x /M± (dashed line) based on Table 5 of Gilfanov etal. 
(2004b). This value is «5-10 times lower than all mean 
values of L x jM- k , suggesting that on average LMXBs play 
a fairly small role in the X-ray emission from our stacked 
samples. Furthermore, late-type galaxies in the local uni- 
verse with similar stellar masses are often found to have 
HMXB emission that is «2-10 times more luminous than 
that expected from LMXBs (see open symbols in Fig. 10b 
and Fig. 3 of Gilfanov et al. 2004c). For galaxy samples se- 
lected via SFR, we find no evidence for significant evolu- 
tion in logLx/SFR (Fig. 10c). We note that if the majority 
of the X-ray flux from our stacked samples originated from 
X-ray-detected sources, then the AGN selection criteria 3 and 
4 from § 4.1 may potentially bias these results. However, we 
find that the more populous and unbiased X-ray-undetected 
source populations dominate the stacked 0.5-8 keV signals at 
the iZ, 80% level (see column 14 of Table 2). For each of the 
three ranges of SFR, logLx/SFR is consistent with a constant 
value, and has a best-fit ratio of logL x /SFR = 39.87 ±0.07 
(X 2 = 13.1 for 12 degrees of freedom). These results sug- 
gest that the contribution from LMXBs is small and that the 
integrated X-ray emission from our late-type galaxies is dom- 
inated by HMXBs. This implies that the evolution of our late- 
type galaxy samples is likely due to changes in star-formation 
activity. 

Since the X-ray emission from our late-type galaxies is 
dominated by star-formation processes, we note that our M+- 
selected stacking results provide a relatively unobscured mea- 
sure of the star-formation activity per unit stellar mass (i.e., 
the SSFR; see Fig. 10b). We find that at z w 1 the X-ray emis- 
sion per unit stellar mass is a factor of f«2-3 larger for galaxies 
with Af* = 3-10 xlO 9 M Q versus that observed for galaxies 
with M* = 1-20 x 10 10 Mq. Atz« 1, we find that L x /M* 
is larger than its local value (S01) by factors of 6 .4 ±2.2 
and 5.8 ±1.6 for late-type galaxies with M* = 3-10 x 10 9 M Q 



T 




log L x (ergs s *) 

FIG. 11. — Logarithm of the star-formation rate log SFR versus the log- 
arithm of the X-ray luminosity logLx f° r normal late-type galaxies. X-ray- 
detected sources from our main sample are indicated as small gray dots. 
Sources that were not detected in the 0.5-8 keV bandpass but were detected 
in either the 0.5-2 keV or 2-8 keV bandpasses are shown as upper limits. 
Results from our X-ray stacking analyses of late-type galaxies selected via 
observed SFR are shown as large filled circles, squares, and triangles, which 
have the same meaning as in Figure 10c; z ~ 3 LBGs that were both un- 
corrected and corrected for AGN contamination have been shown plotted as 
a filled star and open star, respectively. For comparison, we show the local 
galaxy sample from PR07, which includes normal late-type galaxies (open 
circles) and ultraluminous infrared galaxies (ULIRGs; open squares); the 
best-fit PR07 relation is shown as a solid curve. The SFRs for the MW and 
local galaxies M101, M82, and NGC 3256 are indicated (horizontal dashed 
lines). 

and M-t, = 1-20 xl0 10 M Q , respectively. These results are 
broadly consistent with observed differences in the mean SS- 
FRs found by Zheng et al. (2007) for z w 1 galaxies of com- 
parable stellar masses, and imply that the lower-mass galax- 
ies are undergoing more significant stellar mass growth over 
z s=s 0-1 than more massive galaxies. 

In order to quantify the dependences of the X-ray luminos- 
ity on redshift and physical properties, we performed multi- 
variate parametric fitting to our stacked data. For each galaxy 
sample selected via Lb, M+, and SFR, we investigated the red- 
shift evolution of the X-ray luminosity. For this analysis, we 
fit our data to a power-law parametric form: 

log L x (/phys , z) = A + B log /p hys + C log( 1 + z) , (7) 

where / p h ys is a place holder for each of the three physical 
properties (Lb, M+, and SFR) used for our sample selections, 
and A, B, and C are fitting constants. For each sample, we uti- 
lized our X-ray stacking results, equation 7, and x 2 fitting to 
estimate the constants A, B, and C. For our adopted three- 
component model, we constrained A, B, and C using 90% 
confidence errors (Ax 2 = 2.7). The S01 local data points were 
not used for these fits due to differences in galaxy selection, 
instrument calibration, and AGN identification. 

In Table 3, we tabulate our constraints on x 2 , A, B, and C 
for these fits. We find that this particular choice (i.e., eqn. 7) 
of parameterization does not provide acceptable fits for galaxy 
samples selected via L B and M*. However, for galaxy samples 
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FIG. 12. — Cumulative X-ray luminosity dependent AGN fraction 
gt(L2-s kev) versus L 2 -8 keV for late-type galaxies with L B £ 2 X 10 10 Lb,q 
at z 0.6 (filled circles) and z ~ 0.6-1 (open triangles). For reference, we 
have indicated the corresponding AGN fraction for Lyman break galaxies at 
z ~ 3 (filled stars; see § 5 for further details). We have indicated the median 
X-ray detection limit for galaxies in each redshift range (downward-pointing 
arrows along the x-axis). The solid curve represents the best-fit relation for 
gL(Li-% kevX which was fit using all late-type galaxies with z ~ 0-1 and 
L B 2x 1O 1o L B0 in our main sample. For reference, we have shown the 
estimated AGN fraction for galaxies with L B ~ 5 X 10 9 , 10 10 , and 10" L b .q 
(dotted curves; see § 5.2 for further details). 

selected via SFR, we find a good fit for this parameterization 
(X 2 = 9.72 for ten degrees of freedom). We constrain the evo- 
lution of log L x /SFR to be independent of or at most weakly 
dependent on redshift [oc (1 + z )°- 32 ± - 88 ]. 

Based on radio observations of distant star-forming galax- 
ies with SFR w 3-300 M Q yr" 1 in the CDF-N and CDF-S, 
Barger et al. (2007) reported that the X-ray upper limits for 
X-ray-undetected sources were below the level expected from 
the Lx-SFR correlaiton, thus suggesting that the correlation 
may not hold in the high-redshift universe. However, our 
stacking results suggest that the Lx-SFR correlation does 
in fact hold for average galaxies with SFR = 1-5 M yr" 1 
(SFR = 15-100 M Q yr-^outtozwO.S (z~ 1.4). 

For illustrative purposes we have created Figure 11, which 
shows log SFR versus log Lx for normal star-forming galax- 
ies selected from several different sources including X-ray- 
detected galaxies from our main sample (filled gray circles 
and limits), stacked galaxies from this study (filled black 
symbols), local galaxies from PR07 (open circles), local ul- 
traluminous infrared galaxies from PR07 (ULIRGs; open 
squares), and stacked z ~ 3 Lyman break galaxies (LBGs; 
stars; see § 6 for details). For all data used in this plot, we 
have normalized SFRs appropriately to be consistent with our 
adopted Kroupa (2001) IMF and have adjusted X-ray lumi- 
nosities to correspond to the 0.5-8 keV band using a V = 2 
power-law SED. The best-fit Lx-SFR correlations for local 
galaxies from PR07 (solid curve) and z = 0-1.4 late-type 
galaxies from this study (logL x /SFR = 39.87; dotted curve) 
have been shown for reference. For comparison, we have 
shown the SFRs for the MW and local galaxies M101, M82, 
and NGC 3256. 

5.2. AGN Contribution to Stacked Signals 

In this section we estimate the contribution to our stacked 
signals from contaminating AGNs that have X-ray luminosi- 
ties below our X-ray detection limit. This analysis is sim- 
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FIG. 13. — Cumulative redshift-dependent AGN fraction g z (z) ver- 
sus redshift for late-type galaxy samples with L B 10 10 L b ,q and 
^2-8keV ~ 10 41 5 ergs s -1 . We find no significant evolution of g z (z) over 
the redshift range z ~ 0.1-0.8. 

ilar in nature to that in §§ 3.1 and 3.2.2 of Lehmer etal. 
(2007), which was performed for early-type galaxies. We im- 
plement the observed cumulative AGN fraction f c : the frac- 
tion of galaxies harboring an AGN with 2-8 keV luminosity 
of L2-8 kev or greater. Hereafter, we compute f c by taking 
the number of candidate AGNs in a particular galaxy sample 
with a 2-8 keV luminosity of L2-8 kev or greater and dividing 
it by the number of galaxies in which we could have detected 
an AGN with luminosity L2-8 kev- The latter number is com- 
puted by considering the redshift of each galaxy and its corre- 
sponding sensitivity limit, as obtained from spatially varying 
sensitivity maps (see § 4.2 of Alexander et al. 2003 and § 4 of 
Lehmer etal. 2005b); these sensitivity maps were calibrated 
empirically using sources detected by wavdetect at a false- 
positive probability threshold of 1 x 10~ 5 . 

The quantity fc is not only a function of L2-8 kev, but is also 
dependent on the selection of the galaxy sample considered: 
in our case redshift and the physical property / p h ys (i.e., L B , 
M*, and SFR) of the galaxy population may plausibly play a 
role in fc- In our analyses, we assume that each of the respec- 
tive dependencies (i.e., L2-8 kev, z, and / P h ys ) are independent 
of each other, such that f c oc g L (Li-% kev) x g z (z) x gp(/ P h y s), 
where g L , g z , and g p represents the functional dependence 
of the cumulative AGN fraction for each indicated variable 
L2-8 kev, z, and / p h ys , respectively. We made use of the 
2-8 keV bandpass because of its ability to probe relatively 
unattenuated X-ray emission in a regime of the X-ray spec- 
trum where we expect there to be minimal emission from nor- 
mal galaxies (see also criterion 1 of § 4.1 for further details). 
In total 62 («51%) of our 121 X-ray-detected AGN candi- 
dates had 2-8 keV detections; we use these AGNs in our AGN 
fraction analyses. 

We began constructing f c by estimating the shape 
of gL(La-s kev) using late-type galaxies with z = 0-1 and 
Lg £ 2 x 10 10 Lb.q, an optical luminosity regime where we 
have a relatively large number of sources and are sufficiently 
complete out to z = 1 (see Fig Aa). We split this sam- 
ple into two subsets about z = 0.6, to test whether there 
is substantial evolution in the shape and normalization of 
gz.(L2-8kev) over this redshift range. In Figure 12, we 
show gL(Li-% kev) for galaxies in the redshift ranges z £ 0.6 
fcnedian = 0.51; filled circles) and z= 0.6-1 (z me dian = 0.84; 
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FIG. 14. — Cumulative AGN fraction g p (f p h ys ) versus (a)L B , (b)M±, and (c) SFR for late-type galaxies withL 2 _g keV ~ 10 4L5 ergs s -1 at z < 0.6 (dashed curves) 
and z = 0.6-1 (dotted curves). In each plot, we have indicated the best-fit relation (solid line), which is calculated using all galaxies with L2-8 keV ~ 10 4L5 ergs s -1 
and z < 0.6. The downward-pointing arrows indicate which values of / phys were chosen for computing g in equation 8. For reference, the bin sizes of / phys used 
to compute each value of g p (f p hys) and the typical errors of g p (j ph ys ) have been indicated in the lower right-hand corner of each plot. 



open triangles). From Figure 12, we see that the overall shape 
and normalization of gdLi-s kev) for late-type galaxies with 
L B ^2 x 10 Lb,© is similar for galaxies at 

■^median 

= 0.51 

and 

^median — 0.84. We fit the shape of giXLi-ii, kev) using 
least-squares fitting of the 2-8 keV luminosity dependent cu- 
mulative AGN fraction using all galaxies from z = 0-1 with 
L B £ 2 x 10 10 Lg Q. For these fits, we found that the data 
were well-fit by an exponential function, which we parameter- 
ized as logg L (L 2 -8 kev) oc aexpH?(logL 2 -8 keV-39) 2 ], where 
a and b are fitting constants. By construction, this function is 
only valid for logZ^-s kev > 39, which is wl-2 orders of mag- 
nitude less luminous than a typical stacked X-ray luminosity 
of our late-type galaxy samples (see Table 2 and Fig. 10). 
We find best-fit values of a = -0.8 and b = -0.05; in Fig- 
ure 12 {thick black curve), we show our best-fit relation for 
8l(L2-& kev)- 

We constrained further the redshift evolution of fc [i.e., the 
shape of g z (z)] by dividing our main late-type galaxy sample 
into five nearly independent redshift bins (with z = 0.1-0.8) 
and calculating f c for fixed ranges of L2-8 kev and /phys- In 
Figure 13, we show g z (z) as a function of redshift for late-type 
galaxies with L2-8 kev ^ 10 41,5 ergs s" 1 and L B ^ 10 10 L b ,q, 
the approximate completeness limit at z ~ 0.8 (see Fig. Aa). 
Using these data and x 2 fitting [assuming a (1+z)" depen- 
dence], we constrained the redshift evolution of g z (z) to be 
proportional to (1 +z) 015±0 ' 97 ; similar results were found for 
different ranges of L2-8 kev and L B . This result differs from 
that observed for early-type galaxies, where the AGN fraction 
and mean AGN emission has been found to evolve as w(l +z) 3 
(e.g., Brand etal. 2005; Lehmer etal. 2007). 20 We therefore 
conclude that there is little redshift evolution in the late-type 
galaxy AGN fraction over the redshift range z = 0.1-0.8, and 

20 We note that the lack of evolution of the AGN fraction does not nec- 
essarily imply that the optically-luminous (i.e., Lb £ 10 10 Z,s,©) late-type 
galaxy AGN number density is constant over z ~ 0-1. For example, Lg has 
been shown to fade by ssl mag since z as 1 (e.g., Lilly etal. 1995; Wolf etal. 
2003; Faber etal. 2007), suggesting that there are fewer optically luminous 
late-type galaxies in the local universe than at z = 1. 



hereafter we assume that g z (z) remains roughly constant out 
to z = 1 .4. 

To constrain the overall dependence of fc on / p h ys [i.e., 
gpU phys)], we calculated the cumulative AGN fractions for 
late-type galaxy samples with / p h ys by holding the ranges 
of L2-8 kev and z fixed and varying / p h ys - In Figure 14, we 
show gpifphys) versus logL B (Fig. 14a), logM* (Fig. 14b), 
and log SFR (Fig. 14c) for L 2 -8 kev £ 10 41 ' 5 ergs s" 1 at z < 0.6 
(dashed curves) and z= 0.6-1 (dotted curves). We calculated 
gp(fphys) for intervals of / p h ys where we are approximately 
complete at z = 0.6 (for the z ^, 0.6 interval) and z = 1 (for the 
z = 0.6-1 interval; see Fig. 4). For each sample, we again uti- 
lized least-squares fitting to approximate the / p h ys dependence 
°f gpifphy*)- These fits were performed using all data over the 
redshift range z = 0-1 assuming a functional dependence of 
logg^/phys) oc clog/phys- In each panel of Figure 14, we show 
the best-fit solutions for g p (f p hys) with the gray lines. We find 
that the AGN fraction is strongly dependent on L B (Fig. 14a) 
and (Fig. \Ab), such that more optically-luminous and 
massive galaxies have larger AGN fractions; this result is 
consistent with other studies of the AGN host galaxies (e.g., 
Kauffmann etal. 2003b; Nandra etal. 2007; Silverman etal. 
2007). Also, the AGN fraction seems to be mildly dependent 
on the galaxy SFR; however, this is likely due to the fact that 
the SFR is larger on average for more massive galaxies. 

Based on the above estimates of the shapes of gL(L2-& kev), 
g z (z), and gp( fphys), we approximated empirically f c for each 
choice of / p h ys following: 

log/c = logg +aexpH?(logL 2 -8 k eV-39) 2 ]+clog/ phys , (8) 

where go represents the normalization of each relation based 
upon the value of fc at L2-8 kev > 10 41 5 ergs s" 1 and 
L B = 10 10 La.©, M*= 10 10 M Q , and SFR = 2 M Q yr" 1 , for 
galaxy samples selected via L B , M+, and SFR, respectively. 
We note that only go and c are dependent upon our choice 
of / P hys; we find for the set of physical parameters / p h ys = [L B , 
M*, SFR] that g Q = [6.1 x 10" 14 , 8.2 x 10" 8 , 0.20] and c= [1.2, 
0.67, 0.66]. For reference, in Figure 12 we have shown curves 
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of fc for L B = 5x 10 9 , 10 10 , and 10 11 L b .q (dotted curves). 

To estimate the AGN contamination expected for each of 
our stacked samples presented in § 5.1, we followed closely 
the procedure in § 3.2.2 of Lehmer etal. (2007). For com- 
pleteness, we outline this procedure below. 

For each stacked sample, we used equation 8 to compute 
cumulative AGN fractions f c . In Figure 15a, we show f c 
for two of our L B -selected samples: z « 0.2 galaxies with 
Lb = 10 9 5 — 10 10 Lfl.o (dashed curve) and z ~ 1.3 galaxies 
with Lb = 10 10 ' 5 — 10 11 3 Lb.q (solid curve). We then converted 
fc for each sample to a differential AGN fraction (i.e., 
the 2-8 keV luminosity dependent fraction of galaxies har- 
boring AGNs within discrete X-ray luminosity bins of width 
AlogZ.2-8 keV = 0.5; see Fig. 15b). For each sample, we cal- 
culated the 2-8 ke V luminosity dependent fraction of galaxies 
that were below our L2-8 keV detection limit f B (Fig. 15c). We 
then multiplied f D and f B to estimate the 2-8 keV luminosity 
dependent fraction of galaxies harboring an AGN that was un- 
detected in the Chandra exposures fy (i.e., fy = fy x fy; see 
Fig. 15d); these AGNs would not have been removed for our 
stacking analyses. Finally, we approximated the total 2-8 keV 
AGN contamination Lo.5-2 kev(contam) of each stacked sam- 
ple using the following summation: 



L2-& kev(contam) = y~]fu,j X L2-8 keV,i> 



(9) 



where the summation is over all bins of AlogZ.2-8 keV = 0.5 
in the range logZ.2-8 keV = 39^-5. We converted each value of 
L2-& kev(contam) to estimated values of Lo.5-8 kev(contam) and 
Lo.5-2 kev(contam) by assuming the contaminating AGN emis- 
sion roughly follows an X-ray SED described by a power-law. 
In order to constrain the average photon index of the power- 
law, we stacked all X-ray-detected AGNs in our stacked sam- 
ples with 0.5-8 keV luminosities below 10 42 ergs s" 1 . For 
these AGNs, we find a stacked effective photon index of 
Teff = 0.97 ± 0.03. If the intrinsic value of the photon index is 
T = 2, at the median redshift of our main sample (z me dian = 0.5), 
T e ff w 1 corresponds to an intrinsic X-ray column density 
of Na ~ 1-2 xlO 22 cm" 2 . If we assume T e ff = 1 describes 
well the effective SED of the X-ray-undetected AGNs in our 
stacked sample, we find that AGN contamination can account 
for wl-30% (median of w5%) of the 0.5-8 keV emission 
from our stacked samples, suggesting that AGNs are not pro- 
viding a significant contribution to our stacked results. We 
note that the X-ray SED used in this calculation has an impor- 
tant effect on the overall estimate of the AGN contamination. 
Since our estimates for contamination in the 0.5-8 keV and 
0.5-2 keV bandpasses decrease as r e ff decreases, the amount 
of contamination in our samples may be underestimated if our 
choice of T e ff is too flat; however, we find that for conserva- 
tive choices of r e ff (i.e., T e ff £ 2) that are representative of 
even unobscured AGNs (for reference, see Fig. 5), the AGN 
contamination remains low ( 40%) and has no material ef- 
fect on our results. 

It is also worth noting that if there exists a large popula- 
tion of intrinsically luminous yet heavily obscured AGNs in 
our high-redshift galaxy samples that are X-ray-undetected in 
the CDFs, then our conclusions above could be somewhat dif- 
ferent. If such a population were present and had significant 
influence on our stacked results, then we would expect to find 
stacked X-ray spectra that were relatively flat (r e ff £ 1); how- 
ever, as discussed in § 5 . 1 , we find that very few of our stacked 
samples could have such flat spectra in our stacked samples, 
suggesting that such a population, if present, does not have a 
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Fig. 15. — (a) Cumulative AGN fraction (i.e., the fraction of 
galaxies harboring an AGN with a 2-8 keV luminosity of L2-8 kev 
or greater), fc, versus logZ.2-8 kev for two of our stacked samples: 
z ~ 0.2 galaxies with Lb = 10 9 5 — 10 10 Lb.q (dashed curve) and 
z « 1.3 galaxies with L B = 10 10 5 -10 1L3 L B) q (solid curve). Each 
curve was computed following equation 8. (b) Differential AGN 
fractions (i.e., the fraction of galaxies harboring an AGN in discrete 
bins of width AlogZ.2-8 kev = 0.5), fo, versus logZ.2-8 kev- (c) Frac- 
tion of late-type galaxies for which we could not have detected an 
AGN with a 2-8 keV luminosity of L2-8 kev, /b, versus logZ^-s kev- 
(d) Fraction of galaxies harboring AGNs in our optically luminous 
faded samples that would remain undetected due to sensitivity lim- 
itations, fu=f D x f B , versus logZ.2-8 kev; these galaxies would not 
have been removed from our stacking analyses. 

strong effect on our results. 

In Table 2 (cols. 25-27), we have provided the estimated 
fractional AGN contribution to each stacked sample [i.e., 
L £l _ £2 (contam)/L £l _ £2 ] for the FB, SB, and HB using the 
technique described above and an assumed r eff = 1 . We find 
that the estimated AGN contamination is most significant for 
galaxy samples with large values of Lb, M+, and SFR. 

6. EXTENSION TO DISTANT LYMAN BREAK GALAXIES 

As noted in § 1, the global star-formation rate density has 
been observed to increase with redshift out to z ~ 1-1.5. In- 
vestigations of the most-distant LBGs at z ~ 2-7 show that 
the star-formation density peaks around z ~ 1-3 and gradu- 
ally declines toward higher redshifts (e.g., Steidel et al. 1999; 
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FIG. 16. — Logarithm of the X-ray-to-B-band mean luminosity ratio 
log Lx/Lb versus the age of the Universe (for our adopted cosmology, the 
current age of the Universe is 13.47 Gyr) for star-forming galaxies with 
Lg = 3-20 xlO 10 Lb,q. For reference, redshift has been plotted along the 
top axis. We have included mean values of Lx/Lb for the S01 local sam- 
ple of late-type galaxies (open circles), our stacked samples at z ~ 0.5-1.4 
(large filled circles), and z ~ 3 LBGs that were both uncorrected (filled star) 
and corrected (open star) for AGN contamination. At z ~ 0.1-0.4, we have 
plotted Lx/Lb for individual late-type galaxies from our main sample (small 
filled circles and upper limits); the mean values and errors for these galaxies, 
computed using ASURV, has been indicated. 



Giavalisco 2002, 2004b; Bouwens etal. 2004, 2005; Dickin- 
son et al. 2004). To investigate whether the mean X-ray ac- 
tivity from normal late-type galaxies follows a similar trend, 
we study the X-ray properties of a sample of z = 3.01 ± 0.24 
LBGs, which were identified as t/-band "dropouts" through 
the GOODS project (see Lehmer et al. 2005a for details). 

We filtered the GOODS z - 3 LBG sample (Lee et al. 2006) 
to include only LBGs that (1) were within the central «4f0 of 
the «2 Ms CDF-N and «1 Ms CDF-S, and (2) had rest-frame 
B-band luminosities that were similar to the most optically lu- 
minous (L B = 3-20 x 10 10 L b ,q) z = 0-1.4 late-type galaxies 
used in this study, fi-band luminosities were calculated by 
applying fe-corrections to the z 850 -band flux (from GOODS), 
where ^-corrections were derived using an SED appropriate 
for LBGs (see § 2.2 of Lehmer etal. 2005a for details). We 
found that 85 z ~ 3 LBGs from the Lehmer et al. (2005a) sam- 
ple satisfied these two selection criteria. 

We identified three X-ray-detected LBGs, which all had 
L2-8 kev ~ 10 43 ergs s" 1 . Due to their high X-ray luminosi- 
ties, we classified these sources as obvious AGNs. After 
removing these three X-ray-detected AGNs from our z ~ 3 
LBG sample, we performed X-ray stacking analyses as de- 
scribed in § 4.2. We found a significant (3.7a) detection in 
the 0.5-2 keV bandpass, which corresponds roughly to rest- 
frame 2-8 keV emission. Assuming an intrinsic power-law 
X-ray spectrum with a photon index of T = 2, we found a mean 
0.5-8 keV luminosity of L x = (4.9 ± 1.4) x 10 41 ergs s" 1 for 
our z ~ 3 LBGs, a value that agrees well with previous stud- 
ies (e.g., Lehmer et al. 2005a; Laird et al. 2006). 

Since our LBGs reside in the high-redshift universe, AGN 
contamination is expected to have a more significant ef- 
fect on these results than it did for our z = 0-1.4 late-type 
galaxies. The median X-ray luminosity detection limit is 



w8 x 10 41 ergs s" 1 and «10 43 ergs s" 1 for the 0.5-2 keV and 
2-8 keV bandpasses, respectively. To estimate the AGN con- 
tamination, we followed the approach outlined in § 5.2, which 
made use of the 2-8 keV luminosity dependent AGN fraction 



(/c). In Figure 12 (filled stars), we show the X-ray luminos- 
ity dependent AGN fraction gL,LBG(^2-8 kev) for z ^ 3 LBGs. 
We note that at L 2 - 8 kev ~ 10 43 ergs s" 1 , the cumulative AGN 
fraction for z ~ 3 LBGs is a factor of «3-4 times larger than 
that computed for our z = 0-1 late-type galaxy sample with 
similar optical luminosities. Using the functional form for fc 
presented in equation 8, but with a Ri3-4 times larger normal- 
ization factor for z ~ 3 LBGs (i.e., go,LBGs ~ 3-4 go), we find 
that AGN emission may plausibly account for «50-70% of 
the stacked 0.5-2 keV counts. 

In Figure 16, we show the X-ray-to-B-band mean lumi- 
nosity ratio for L B = 3-20 x 10 10 L B , Q star-forming galaxies 
(i.e., late-type galaxies and LBGs) as a function of the age of 
the Universe. Together, these data span «85% of cosmic his- 
tory (i.e., out to z ~ 3). As presented in § 5.1, L X IL B shows 
significant evolution over the redshift range z = 0-1 .4, and af- 
ter correcting for AGN contamination, we find that L X /L B is 
similar for z ~ 3 LBGs and z = 1 .4 late-type galaxies. This 
result suggests that the non-AGN X-ray emission for the most 
luminous star-forming galaxies may flatten near z ~ 1 .4-3, 
which has been predicted roughly from simulations of how the 
normal-galaxy X-ray emission is expected to respond due to 
global changes in the star-formation rate density (e.g., Ghosh 
& White 2001). 

To test whether the Lx-SFR correlation is similar for 
z ~ 3 LBGs as we found for late-type galaxies at z = 0-1.4, 
we approximated absorption-corrected SFRs for the z <~ 3 
LBGs using UV band emission. These SFRs were ap- 
proximated following SFR = 9.8 x 10~ n 7L uv , where 7 « 6 
is the absorption-correction factor (see Giavalisco etal. 
2004b). As described in § 3.3, we approximated the UV 
luminosity following Ltjv = 3.3fZ„(2800 A); however, here 
/„(2800 A) was derived using our adopted LBG SED. We 
find that the mean SFR for optically-luminous z ~ 3 LBGs 
is «60M Q yr" 1 . After correcting the mean stacked X-ray 
luminosity for AGN contamination, we find an X-ray-to- 
SFR ratio of logL x /SFR w 39.4-39.6, which is suggestively 
lower than its value for z= 0-1.4 late-type galaxies (i.e., 
logLx/SFR = 39.87 ± 0.04). We note, however, that the esti- 
mation of the mean SFR for these galaxies is based solely on 
the 2800 A emission and is therefore highly uncertain; further 
constraints on the mean infrared luminosity of these sources 
would help considerably. For reference, in Figure 1 1 we have 
plotted log SFR versus logLx for the z ~ 3 LBGs both uncor- 
rected (filled star) and corrected (open star) for AGN contam- 
ination. 

It is interesting to note that the AGN-corrected LBG sam- 
ple appears to have L x and SFR values similar to those of 
local ULIRGs, which are relatively X-ray underluminous for 
their derived SFRs. PR07 suggest that galaxies with SFRs that 
begin to exceed >, 50-100 M Q yr" 1 , similar to the ULIRGs, 
may have X-ray properties that are completely dominated by 
HMXBs. In such systems, other X-ray-emitting populations 
that may normally be significant in more quiescent galaxies 
(e.g., LMXBs, hot gas, supernovae and their remnants, etc.) 
would collectively contribute only a negligible fraction of the 
total X-ray emission, thus causing the overall galaxies to ap- 
pear X-ray underluminous compared with the Lx-SFR corre- 
lation. We note, however, that our results presented in § 5.1 
suggest that HMXBs likely dominate the X-ray emission from 
galaxies with SFRs much lower than those of ULIRGs. An- 
other possibility is that as the absorption within the most ac- 
tively star-forming galaxies increases with SFR (see discus- 
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sion in § 5.1), the X-ray emission from point-source popula- 
tions may become significantly obscured. Deep X-ray obser- 
vations of local LIRGs and ULIRGs with Chandra and XMM- 
Newton have found significant X-ray absorption in their point- 
source populations (see, e.g., Lira etal. 2002; Zezas etal. 
2002, 2003) thus providing some support for this possibility. 
Future investigations into the nature of the X-ray populations 
of ULIRGs could help resolve this issue. 

7. SUMMARY AND FUTURE WORK 

We have investigated the X-ray emission from 2568 nor- 
mal late-type galaxies over the redshift range z = 0-1 .4 that lie 
within the Chandra Deep Fields (CDFs). Our late-type galaxy 
sample was constructed primarily using color-magnitude dia- 
grams, which incorporated rest-frame U — V color and abso- 
lute V-band magnitudes, to isolate blue late-type galaxies (see 
§ 2 for details). In total, 225 («9%) of our late-type galax- 
ies were detected individually in the X-ray band. Based on 
X-ray and optical spectral properties, X-ray-to-optical flux 
ratios, the correlation between X-ray luminosity and star- 
formation rate, and comparisons between infrared and radio 
properties, we infer that 121 («53%) of the X-ray-detected 
late-type galaxies are dominated by AGN emission. The re- 
maining 104 X-ray-detected galaxies had X-ray and multi- 
wavelength properties consistent with those of normal late- 
type galaxies with X-ray emission dominated by X-ray bina- 
ries (HMXBs and LMXBs). To study the X-ray emission and 
evolution from large representative populations of late-type 
galaxies (i.e., including the X-ray-undetected sources), we 
utilized X-ray stacking analyses of galaxy populations with 
AGN candidates removed. We stacked normal-galaxy sam- 
ples that were selected via their rest-frame fi-band luminosity 
(Lb), stellar mass (M*), and star-formation rate (SFR) in red- 
shift bins (see § 5.1). Furthermore, we compared these results 
with those found for a sample of z ~ 3 LBGs from Lehmer 
etal. (2005a). In the points below, we summarize our key 
findings: 

1 . We obtained significant detections in the 0.5-2 keV and 
0.5-8 keV bandpasses for all of our stacked samples. We es- 
timated that LMXBs and low-level AGNs provide only low- 
level contributions to the stacked X-ray emission from our 
samples and that HMXBs constitute the dominant X-ray- 
emitting component. Therefore, for these galaxies the X-ray 
emission is tracing primarily star-formation activity. Nor- 
mal late-type galaxy samples selected via Lb and M t show 
significant (at the >99.9% confidence level) evolution in 
their average X-ray properties from z = to 1.4. For the 
most optically luminous (Lb ~ 5 x 10 10 Lb.q) and massive 
(M* «3x 10 10 Mq) late-type galaxies at z = 1 .4, L X /L B and 
Lx/M* are measured to be larger than their local values by 
factors of ^2-4 and w7-13, respectively. 

2. We find that late-type galaxies of lower stellar mass gen- 
erally have larger X-ray-to-stellar-mass mean ratios (Lx/M+) 
than their higher-mass analogs. Over z ~ 0.2-1, galaxies 
with M 4 «6x 10 9 M Q are a factor of w2-3 times more 
X-ray luminous per unit stellar mass than galaxies with M+ 

3 x 10 10 Mq. 

3. We characterized the X-ray properties of 888 24/im- 
detected late-type galaxies in our sample. The 24/im data 
allowed us to select galaxy samples selected via their SFRs. 
For these samples, we found that the X-ray luminosity is well- 
predicted by a constant Lx-to-SFR ratio, similar to the Lx- 
SFR correlation reported by previous authors (e.g., PR07). 



This implies that the Lx-SFR correlation holds out to at least 
z = 0.5, 1, and 1.4 for galaxies with SFR » 2, 10, and 50 M Q 
yr" 1 , respectively, and supports the idea that the strong X-ray 
evolution observed for normal late-type galaxies selected via 
Lb and is likely due to strong changes in SFR. 

4. The X-ray properties of our most optically-luminous 
(Lb = 3-20 xl0*° L#q) late-type galaxies at z = 1.4 are 
comparable to those for z ~ 3 LBGs with similar optical 
luminosities, once X-ray-undetected AGN contamination in 
the LBG population has been accounted for. This suggests 
that there may plausibly be a flattening in the Lx/Lb-z 
relation for optically-luminous star-forming galaxies between 
Z ~ 1 .4-3. We estimate a mean SFR of «60 Mq yr" 1 for these 
LBGs. We find that the observed mean X-ray luminosity is 
suggestively underluminous based on the Lx-SFR correlation 
prediction; this result is similar to that found for local 
ULIRGs with comparable SFRs. 

The above results can be improved greatly through (1) the 
study of late-type galaxy populations in other existing multi- 
wavelength extragalactic survey fields that contain Chandra 
observations, (2) additional observations of the CDFs to even 
better sensitivity levels than are currently available, and (3) 
observations with future X-ray missions with imaging capa- 
bilities that are complementary to those of Chandra. 

The areal footprint of the CDF regions used in this study 
totals «0.18 deg (see § 2.1 for details), which severely lim- 
its the number of galaxies being studied at z ,$ 0.2. The red- 
shift range z = 0-0.2 spans w2.5 Gyr of cosmic look-back 
time, compared with the «6.6 Gyr over z = 0.2-1.4, that we 
are most effectively studying above. Several complemen- 
tary deep and wide extragalactic Chandra survey fields that 
have recently been conducted or are in progress can improve 
the present situation. Three such ideal survey fields are the 
w200 ks All-wavelength Extended Groth Strip International 
Survey (AEGIS; e.g., Nandra etal. 2005; Davis etal. 2006), 
the «50 ks Chandra Cosmic Evolution Survey (C-COSMOS; 
PI: M. Elvis; see also, Scoville etal. 2007), 21 and the «5 ks 
NOAO Deep Wide-Field Survey (NDWFS; e.g., Murray et al. 
2005). 

Our study shows that identifying and removing contaminat- 
ing AGNs is a crucial ingredient to studying the stacked X-ray 
properties of normal-galaxy samples (see § 5.2). In order to 
avoid significant levels of AGN contamination in shallower 
wide-area Chandra surveys such as AEGIS, C-COSMOS, and 
NDWFS, X-ray studies of normal galaxy populations must 
be limited to lower redshift intervals. If we require a de- 
tection limit of Lo.5-2 keV = 10 41 ' 5 ergs s _I over regions of 
each survey where the Chandra sensitivity is optimal (i.e., 
6.'0 from the Chandra aim-points), then we estimate that 
AEGIS, C-COSMOS, and NDWFS could effectively be used 
to study normal late-type galaxy populations at redshifts less 
than «0.5, «0.25, and ssO.l, respectively. When factoring in 
the areal coverage of these fields, we estimate that AEGIS and 
C-COSMOS will contain w250-350 late-type galaxies with 

21 We note that the C-COSMOS observing strategy has been designed to 
have several overlapping f»50 ks ACIS-I exposures; in these regions, the to- 
tal exposure will reach fs:200 ks due to the overlaping exposures. However, 
the analyses presented in this paper utilize only high-quality Chandra imag- 
ing at off-axis angles of 6'0 (see § 2.1 for additional justification), where 
there is very little overlap in the C-COSMOS Chandra pointings. We there- 
fore do not consider the overlapping regions of the C-COSMOS exposures 
in this discussion. For additional information regarding C-COSMOS, see 
http://cfa-www.harvard.edu/hea/cos/C-COSMOS.html 
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L B > 3 x 10 9 L B q at z = 0.1-0.2 (i.e., a factor of «3 larger 
than the CDFs) and that AEGIS, C-COSMOS, and NDWFS 
taken together will have w 100-250 similar galaxies at z ^ 0.1 
(i.e., a factor of sa22 times larger than the CDFs). Study- 
ing these galaxy populations would significantly improve con- 
straints on the X-ray properties of z ~ 0.2 late-type galaxies. 

Deeper Chandra observations of the CDFs would provide 
additional insight into the X-ray properties of the normal 
late-type galaxy populations presented in this study and en- 
able us to extend our analyses to higher redshifts. Typi- 
cal galaxies in our sample have mean X-ray luminosities of 
L x ~ 10 39 5 -10 40 ergs s" 1 . In the «2 Ms CDF-N, where our 
Chandra sensitivity is greatest, we expect that galaxies with 
Lx ~ 10 395 ergs s" 1 should be detectable out to z ~ 0.2. Out 
of the 1 1 late-type galaxies with Lb £ 10 9 Lb.q atz £ 0.2, we 
find X-ray detections for 8 of them (w72%). Deeper Chan- 
dra observations over the CDFs will allow for the detection 
of such galaxies out to progressively higher redshifts. For 
Chandra exposures of «5 Ms and wlO Ms, a source with 
Lx ~ 10 39 5 ergs s" 1 should be detectable out to z ~ 0.3 and 
z ~ 0.4, respectively. If we pessimistically assume an X-ray 
detection fraction of sa70%, we estimate that «30 and «50 
normal galaxies would be detected individually for each re- 
spective deep exposure. The individual detection of these 
galaxies would allow for improved constraints on both the 
X-ray luminosity functions of normal late-type galaxies at 
higher redshifts and the low-luminosity AGN contributions 
to our stacked signals. At higher redshifts, an «5 Ms and 
wlO Ms Chandra exposure would enable us to effectively 
perform studies of normal late-type galaxies, similar to those 
presented in this paper, out to z ~ 2 and z ~ 2.5, respectively. 
Such a data set would provide, for the first time, a reliable 
constraint on the X-ray emission from normal star-forming 
galaxies near the peak of the global star-formation rate den- 

22 For additional information regarding the future X-ray mis- 
sions XEUS and Generation-X, see http://www.rssd.esa.int/XEUS/ and 
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sity at z ~ 1.5-3. 

In addition to the improvement that additional Chandra 
observations could provide, future X-ray missions such as 
XEUS and Generation-X 22 should enable new scientific in- 
vestigations of distant normal galaxies. XEUS will be able 
to place significant spectral constraints for sources with 
0.5-2 keV fluxes down to «10~ 17 ergs cm" 2 s" 1 , a level fainter 
than many of the average fluxes derived from our samples. 
Generation-X is planned to provide imaging with «0."1 per 
resolution element and will easily probe to flux levels of 
wlO" 19 ergs cm" 2 s" 1 (L x ~ 5 x 10 38 ergs s" 1 at z = 1). At 
these levels, the first detailed investigations of the evolution 
of the normal late-type galaxy X-ray luminosity function can 
be performed effectively out to z 1 without confusion prob- 
lems due to the crowding of large numbers of sources. Fur- 
thermore, these observations will allow for new constraints to 
be placed on the X-ray populations of all normal galaxies in 
the observable universe that are offset by more than w0.8 kpc 
from their host-galaxy nuclei. 
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